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1.  INTROPUCnON 


1.1  BACKGROUND 

Frozen  soils  exist  in  Arctic  regions  as  perennially  frozen  soil  (permafrost)  that 
underlies  nearly  20%  of  the  land  surface  of  the  earth.  In  addition,  artificial  freezing  is 
widely  used  as  a  construction  aid  to  enhance  ground  support  and  for  control  of  water 
during  excavations  of  tunnels,  shafts,  etc.  Projects  involving  construction  in  or  on  frozen 
earth  pose  unique  design  problems  because  the  time-dependent  strength-deformation 
behavior  of  frozen  soil  is  probably  the  most  variable  and  difficult  of  all  "geomaterials"  to 
understand  and  model.  Even  the  simplest  system,  such  as  a  saturated  frozen  sand, 
represents  a  highly  complex  interaction  between  the  solid  soil  skeleton  and  the  pore  matrix, 
composed  of  ice  and  unfrozen  water,  that  changes  continuously  with  time  as  a  function  of 
temperature,  state  of  stress  and/or  strain  level. 

In  order  to  develop  solutions  to  design  problems  involving  frozen  soils,  engineers 
have  generally  relied  on  empirically  based  behavioral  models  as  a  means  of  predicting 
material  response.  These  have  been  based  on  conelations  between  various  parameters  such 
as  temperature,  strain  rate  and  peak  strength,  or  creep  stress  and  minimum  creep  rate  or 
time  to  minimum  creep  rate.  Unconfined  compression  tests  have  been  used  most  often  as  a 
means  of  calibrating  these  behavioral  models.  These  tests  are  performed  in  the  laboratory 
where  the  time  scale  is  generally  much  shorter  than  the  field.  Accordingly,  extrapolation  of 
the  measured  laboratory  behavior  to  the  field  can  introduce  a  considerable  amount  of 
uncertainty  in  the  predicted  response. 

Unconfined  compression  tests  also  apply  a  relatively  simple  state  of  stress  to  the 
specimen.  In  the  field,  where  the  actual  states  of  stress  and  deformational  constraints  are 
very  complex,  behavior  extrapolated  from  such  simple  element  tests  may  lead  to 
unconservative  predictions  unless  a  large  factor  of  safety  (factor  of  ignorance)  is  used.  In 
order  to  model  complex  field  conditions,  it  would  be  necessary  to  conduct  sophisticated 
laboratory  testing  programs  where  principal  stress  directions  could  be  controlled  and 
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widdy  varying  deformational  constraints  could  be  imposed.  Such  testing  procedures  are 
extremely  costly  and  have  not  been  developed  for  frozen  soil  applications. 

Most  experimental  data  on  frozen  soils  comes  from  laboratory  tests  wherein  the 
specimens  have  been  prepared  vrith  "zero”  effective  (intergranular)  stress  acting  on  the  soil 
skdeton  prior  to  freezing.  In  triaxial  compression  tests,  a  confining  pressure  is 
subsequently  applied  to  the  frozen  soil  and,  if  it  is  completdy  saturated,  this  confinement 
is  initially  carried  primarily  by  the  pore  ice.  In  contrast,  both  natural  and  artificial 
freezing  occurs  in  soils  that  already  have  significant  effective  stresses  acting  on  the  soil 
skdeton  prior  to  freezing.  No  data  exist  to  predict  how  differences  in  the  pre— freezing 
effective  stresses  may  affect  laboratory  versus  fidd  behavior. 

Therefore,  empirically  based  behavioral  models,  while  certainly  necessary  for  solving 
immediate  engineering  design  problems,  are  liable  to  involve  significant  uncertainty  in  thdr 
ability  to  predict  actual  behavior  in  the  fidd.  The  eventual  devdopment  of  more 
rationally  based  design  procedures,  which  presumably  can  lead  to  safer  and  more  cost 
efficient  designs,  requires  a  far  better  understanding  of  the  physical  mechanisms  that 
control  the  time-dependent  behavior  of  frozen  soils. 

1.2  OBJECTIVES  OP  RESEARCH 

The  long  term  objective  of  MIT  frozen  soil  research  is  to  devdop  physically  based 
constitutive  relations  for  frozen  soil  behavior,  for  which  an  understanding  of  the 
fundamental  behavior  of  frozen  soils  is  essential.  The  near  term  research  objective  is  the 
identification  of  the  principal  mechanisms  which  control  the  behavior  of  frozen  soils  and  a 
quantitative  assessment  of  their  relative  importance. 

These  mechanisms  can  be  broadly  classified  into  those  which  are  associated  with  the 
soil  skdeton,  those  associated  with  the  pore  ice,  and  those  which  are  due  to  the  interaction 
between  the  soil  partides  and  the  pore  ice.  Ideally  frozen  soils  could  be  studied  on  the 
micro-level.  Direct  measurements  of  the  behavior  of  the  soil  skeleton  independent  of  the 
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ice  matrix  and  visa  versa  would  go  a  long  way  towards  understanding  the  complex 
interaction  between  these  two  components.  But  such  measurements  are  not  currently 
possible.  An  alternative  approach  is  to  study  frozen  soils  at  a  macro-level.  This  involves 
observing  the  overall  behavior  of  the  frozen  soil  system  under  a  wide  variety  of  conditions 
and  comparing  it  to  predicted  behavior  using  simplified  models  aimed  at  describing  the 
various  physical  mechanisms.  The  relative  importance  of  these  mechanisms  can  then  be 
addressed.  The  primary  trust  of  MIT’s  current  research  has  taken  this  macro-level 
approach. 

1.2.1  Rational  and  Original  Objectives 

Ting  et  al.  (1983)  proposed  and  attempted  to  evaluate  various  physical 
mechanisms  controlling  the  strength  of  frozen  sand,  concluding  that  the  principal 
components  are  (see  Table  1.1):  ice  strength  and  strengthening;  soil  strength;  and  soil 
skeleton-ice  matrix  interaction  through  dilatancy  effects  and  structural  hindrance.  Figure 
1.1(a)  shows  their  proposed  structure  for  frozen  sand  and  Fig.  1.1(b)  illustrates  their 
attempt  to  quantify  the  relative  importance  of  these  different  strength  mechanisms  as  a 
function  of  the  volume  fraction  of  sand  from  unconfined  compression  data  obtained  by 
Goughnour  and  Andersland  (1968)  on  frozen  Ottawa  sand.  A  major  objective  of  our  1988 
ARO  proposal  was  to  extend  the  Ting  et  al.  conceptual  picture  to  obtain  a  quantitative 
assessment  of  the  physical  components  of  strength  for  a  range  of  frozen  soils  and  testing 
conditions.  Specifically,  the  research  would  make,  for  the  first  time,  a  direct  comparison  of 
the  total  stress  versus  strain  behavior  of  representative  frozen  soils  with  the  undrained 
total  and  effective  stress  versus  strain  behavior  of  the  same  systems  in  an  unfrozen  state. 
By  varying  the  testing  conditions  (i.e.,  specimen  density,  confining  stress,  strain  rate,  etc.) 
for  each  system,  the  research  hoped  to  quantify  the  factors  affecting  the  various  physical 
mechanisms  and  their  relative  importance.  The  results  would  also  be  compared  to 
predictions  by  Ladanyi’s  "dilatancy— hardening”  theory  (e.g.,  Ladanyi  and  Morel  1990). 
This  theory  is  further  described  in  Section  5  and  assumes  that  the  strength  of  frozen  dense 
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sand  has  two  components:  1)  one  due  to  the  undrained  strength  of  the  soil  skeleton,  which 
is  controlled  by  its  friction  angle  and  the  tensile  strength  of  ice  (at  low  confining  stresses); 
and  2)  the  other  due  to  the  compressive  strength  of  the  ice  matrix  (equated  to  that  of 
regular  granular  ice  at  the  same  temperature,  strain  rate  and  confining  pressure). 

Limited,  but  consistent,  volumetric  strain  data  in  the  literature  indicated  that 
shear  of  frozen  soils  in  triaxial  compression  occurred  at  almost  constant  volume  up  to  the 
peak  strength.  This  formed  the  basis  for  Ladanyi’s  hypothesis  that  the  undrained  (AV  = 
0)  effective  stress-strain  behavior  of  the  same  soil  in  an  unfrozen  state  could  provide  an 
estimate  of  the  contribution  of  the  soil  skeleton  strength  (components  2  and  3bi  in  Table 
1.1)  to  the  overall  strength  of  the  frozen  soil.  Consequently,  the  experimental  program 
would  conduct  isotropically  consolidated— undrained  triaxial  compression  (CIUC)  tests  with 
pore  pressure  measurements  on  unfrozen  soil  in  order  to  compare  its  effective  stress  and 
total  stress  behavior  with  the  total  stress  behavior  of  the  same  soil  in  a  frozen  state.  Table 
1.2  shows  the  proposed  scope  of  a  three  year  strength  testing  program.  Some  explanations 
follow: 

1)  The  program  would  start  with  Manchester  Fine  Sand  since  its  behavior  is  typical  of 
frozen  sands  and  prior  ARO  sponsored  research  (Martin  et  al.  1981;  Ting  et  al. 
1983)  had  developed  extensive  unconfined  creep  data  at  varying  densities  and 
degrees  of  ice  saturation. 

2)  The  confining  stresses  for  all  test  series  would  vary  over  a  vdde  range  (0.1  to  10 
MPa),  since  this  is  an  important  variable  for  both  unfrozen  and  frozen  soil.  All 
tests  would  be  conducted  with  measurements  of  volumetric  strain  and  would 
incorporate  lubricated  end  caps. 

3)  The  U  Series  would  represent  conventional  CIUC  (isotropically  consolidated— 
undrained  compression)  tests  on  unfrozen  soil  systems. 

4)  The  F  Series  would  represent  conventional  UUC  tests  on  frozen  soil  (i.e.,  specimens 
are  sheared  immediately  after  application  of  the  cell  pressure,  ac);  the  FC  Series 
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would  involye  very  prolonged  application  of  confining  stress  in  order  to  (x>nsolIdate 
the  frozen  soil  (i.e.,  allow  ice  to  flow  out  of  the  specimen  until  the  soil  skeleton 
carries  the  confining  stress). 

5)  The  FI  Series  would  evaluate  the  influence  of  confining  stress  at  a  typical 
temperature  and  strain  rate.  The  F2  Series  would  evaluate  the  influence  of  higher 
and  lower  values  of  temperature  and  e  at  low  and  high  confining  pressures. 

6)  The  frozen  Boston  Blue  Clay  would  represent  a  soil  having  a  strength 
approximately  equal  to  that  of  granular  ice  and  frozen  Fairbanks  Silt  would 
represent  a  soil  having  a  strength  significantly  less  than  granular  ice. 

7)  The  program  on  Glass  Beads  was  included  for  comparison  with  Ting  et  al.  (1983) 
aeep  data. 

1.2.2  Actual  Scope  and  Objectives 

The  actual  scope  of  the  triaxial  testing  progra.m  has  been  restricted  to  the  U 
and  F  Series  on  Manchester  Fine  Sand  (MFS)  and  hence  has  far  less  breadth  than 
originally  proposed.  Some  explanations  are  needed  to  understand  the  reasons  for  this 
substantial  change. 

Development  of  New  Experimental  Capabilities 
The  time  and  effort  required  to  design,  construct  and  evaluate  two  new  triaxial 
systems  (especially  that  for  frozen  soil)  were  greater  than  expected  in  large  part  due  to 
changes  in  measurement  and  experimental  techniques  as  follows: 

•  The  addition  of  on— specimen  measurements  of  axial  strain  (which  required  a  new 
design  since  suitable  devices  were  not  commercially  available).  This  proved  to  be 
both  essential  and  highly  successful. 

•  The  subsequent  need  to  develop  a  complex  technique  to  measure  volumetric  strains 
via  changes  in  cell  volume  rather  than  radial  strain  gages  as  originally  proposed. 
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•  The  change  from  "wet  tamping"  to  "multiple  sieve  pluviation"  to  achieve  more 
uniform  specimen  densities. 

Consequently  it  took  about  one  year  of  full  time  effort  by  one  student  (Glen  Andersen)  to 
produce  the  first  "reported"  test  on  frozen  sand.  A  second  student  (Chris  Swan)  devoted  a 
similar  effort  for  the  unfrozen  sand  testing  system. 

Results  from  Testing  Frozen  MFS 

The  improved  measurement  techniques  led  to  the  acquisition  of  unique  stress— 
strain-strength  data  compared  to  that  reported  in  the  literature,  particularly  regarding 
behavior  at  small  strains.  The  addition  of  on-spedmen  axial  strain  devices  allowed  for  the 
measurement  of  strains  down  to  less  than  0.01%.  Consequently,  the  decision  was  made  to 
obtain  for  the  Rrst  time  the  complete  characterization  of  the  stress-strain  behavior  of  a 
frozen  sand  in  triaxial  compression  as  a  function  of  relative  density  (Dr),  confining  stress 
(ffc),  strain  rate  («)  and  temperature  (T). 

Research  Objectives 

The  main  objective  of  the  research  was  to  use  the  new  experimental  system  in 
order  to  obtain  detailed  measurements  of  the  behavior  of  frozen  Manchester  Fine  Sand 
(MFS)  as  a  function  of  the  above  four  variables  (i.e.,  density,  confinement,  strain  rate  and 
temperature).  An  extensive  literature  review  revealed  that  information  of  this  type  had 
not  been  published  for  any  frozen  soil.  The  authors  feel  that  this  complete  characterization 
represents  the  most  significant  contribution  of  the  research  and  is  a  necessary  and 
important  step  towards  the  eventual  development  of  physically  based  constitutive  models 
for  frozen  soils. 

The  second  objective  of  the  research  was  to  evaluate  the  measured  behavior  of 
frozen  Manchester  Fine  Sand  (MFS)  in  terms  of  Ladanyi’s  (1985)  dilatancy-hardening 
model.  His  model  attempts  to  quantify  the  relative  contributions  of  the  frictional 
resistance  of  the  soil  skeleton  and  of  the  tensile  and  compressive  strength  of  the  ice  matrix. 
The  evaluation  process  required  detailed  measurements  of  the  stress-strain-strength 
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behavior  of  un&ozen  MFS  over  the  same  range  of  densities  and  confining  pressures  as  used 
for  the  test  series  on  frozen  sand.  This  effort  also  constitutes  a  significant  contribution  and 
the  first  data  set  comparing  undrained— unfrozen  sand  and  frozen  sand  behavior  over  a  wide 
range  of  densities,  strain  rates  and  confining  pressures. 

The  third  and  final  objective  evolved  during  the  research  and  involved  the  use 
of  a  behavioral  model  from  the  field  of  composite  materials  in  an  effort  to  understand  the 
interaction  between  the  sand  skeleton  and  pore  ice  at  small  strains.  In  particular, 
predictions  from  Gounto’s  (1964)  composite  material  model  for  concrete  are  compared  to 
measured  behavior. 

1.3  ORGANIZATION  OF  REPORT 

Section  2  describes  the  laboratory  facilities  and  equipment  and  triaxial  testing 
procedures  used  in  both  frozen  and  unfrozen  tests.  It  also  contains  a  description  of 
Manchester  Fine  Sand  and  specimen  preparation  techniques. 

Section  3  presents  the  measured  behavior  of  frozen  Manchester  Fine  Sand  in  triaxial 
compression.  It  gives  an  overview  of  the  scope  of  the  testing  program  and  then  shows  the 
effects  of  the  relative  density,  confining  pressure  and  strain  rate  on  the  stress-strain 
behavior  and  is  based  on  data  contained  in  the  doctoral  thesis  by  Andersen  (1991).  These 
tests  were  conducted  at  one  temperature  (-9.5*  C).  Summary  plots  of  parameters  taken 
directly  from  the  stress-strain  curves  are  presented  and  form  the  basis  of  discussion.  These 
plots  are  supplemented  by  typical  stress  and  volumetric  strain  versus  axial  strain  curves. 
Parameters  are  tabulated  in  Appendix  A.  The  effect  of  temperature,  which  is  still  being 
studied,  and  the  influence  of  different  specimen  preparation  techniques  on  the  behavior  of 
frozen  MFS  are  also  presented.  In  addition,  the  measured  behavior  at  very  small  strains  is 
analyzed  in  terms  of  the  composite  materials  model  presented  by  Counto  (1964). 

The  results  of  triaxial  tests  on  unfrozen  MFS  are  presented  in  Section  4.  After 
outlining  the  scope  of  the  testing,  the  effects  of  relative  density  and  effective  confining 
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pressure  on  nndrained  stress-strain  behavior  are  presented.  Tabulated  test  results  are 
presented  in  Appendix  B.  The  unfrozen  test  results  are  analyzed  using  the  "state 
parameter"  concept  developed  by  Been  and  Jefferies  (1985). 

Section  5  presents  the  dilantanc^— hardening  modd  devdoped  by  Ladanyi  (1985). 
Predicted  results  from  this  modd  are  compared  to  those  measured  in  this  research. 

Section  6  summarizes  the  prindpal  condusions  and  presents  recommendations  for 
continued  research. 

1.4  PUBLICATIONS  AND  REPORTS 

Progress  Reports  were  submitted  to  ARO  for  the  six  month  periods  ending  in  June 
and  December  of  1989,  1990  and  1991.  A  copy  of  the  doctoral  thesis  that  was  completed  in 
June,  1991  by  Glen  R.  Andersen  entitled  "Physical  Mechanisms  Controlling  the  Strength 
and  Deformation  Behavior  of  Frozen  Sand"  was  sent  to  ARO.  Dr.  Andersen  is  currently 
preparing  drafts  of  two  co-authored  papers  for  submittal  to  the  Canadian  Geotechnical 
Journal.  The  first  will  present  the  effects  of  density,  confining  pressure  and  strain  rate 
(and  maybe  temperature)  on  the  behavior  of  ffozen  MFS.  The  second  will  analyze  the  data 
(e.g.,  use  of  composite  and  Ladanyi’s  models)  and  compare  these  data  with  prior  results. 
(Note:  the  papers  will  be  too  long  for  acceptance  by  ASCE,  whereas  the  CGJ  has  no  page 
limitations  and  has  contained  many  excellent  papas  on  ice  and  frozen  soils). 

1.5  STAFF 

Dr.  Charles  C.  Ladd,  Professor  of  Civil  Engineering,  and  Dr.  John  T.  Germaine, 
Principal  Research  Associate,  served  as  co-principal  investigators  for  the  three  year 
project.  In  simplified  terms.  Dr.  Germaine  was  mainly  responsible  for  development  of  the 
new  testing  capabilities  and  supervision  of  the  experimental  work,  whereas  Professor  Ladd 
administered  the  contract  and  was  mainly  responsible  for  supervision  of  the  project  and  the 
writing  of  reports  and  the  thesis. 


8 


Glen  R.  Andersen  worked  on  the  project  for  two  and  one-half  years  either  as  a 
full— time  graduate  Research  Assistant  or  a  part-time  Instructor,  made  an  extensive  review 
of  the  literature  on  the  behavior  of  i(%  and  frozen  sands,  and  conducted  or  supervised  most 
of  the  tests  on  frozen  MFS  involving  the  effects  of  density,  confining  pressure  and  strain 
rate.  After  receiving  his  ScD  degree  in  June  1991,  Dr.  Andersen  joined  the  Chevron  Oil 
Field  Research  Company  in  La  Habra,  California. 

Chris  W.  Swan  worked  on  the  project  through  the  last  two  years,  first  as  a  recipient 
of  a  Department  of  Civil  Engineering  Fellowship  and  then  as  a  full-time  graduate  Research 
Assistant.  He  developed  the  capability  for  conducting  high  pressure  triaxial  tests  on 
unfrozen  MFS  that  produced  the  results  presented  in  Section  4  and  then  assumed 
responsibility  for  testing  frozen  MFS,  primarily  on  the  effects  of  temperature.  Mr.  Swan 
plans  to  continue  research  in  this  area  for  his  doctoral  thesis,  with  emphasis  on  comparison 
of  the  influence  of  the  magnitude  of  the  pre-freezing  effective  stress  on  the  measured 
behavior  of  frozen  MFS. 

Other  involvement  by  students  included: 

1)  Patrick  Kinnicut  as  a  half-time  graduate  Research  Assistant  during  the  Spring 
semester  of  1989  who  developed  improved  techniques  for  preparation  of  specimens  of 
MFS; 

2)  Peter  Kukielski  as  a  Summer  Intern  Fellowship  Student  from  the  University  of 
Rhode  Island  who  performed  a  significant  number  of  high  quality  tests  on  frozen 
MFS  under  Mr.  Andersen’s  supervision; 

3)  ^ctor  Salvador  as  a  Summer  Intern  Fellowship  Student  from  Columbia  University 
who  assisted  Mr.  Swan  in  conducting  tests  on  frozen  MFS  at  varying  temperatures. 
Mr.  David  M.  Cole  of  the  Cold  Regions  Research  and  Engineering  Laboratory 

(CRUEL)  in  Hanover,  NH  served  as  technical  liaison  for  the  project.  He  also  served  as  a 
member  of  Andersen’s  doctoral  thesis  committee  along  with  Professor  Ladd  and  Dr. 
Germaine. 
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(c)  Tension  in  unfrozen  water  film  •  Probably  not  very  significant  (e.g.,  decreasing  ice 

saturation  leads  to  much  lower  strength). 


TabteO  TettUtiTeSa^ecfTdttdalStieiigtliTeatiiig 
Piognm  ia  ARO  Pxoposal 


A.  TEST  SERIES  (Stresses  in  MPa  and  e  per  second) 


Designation 

and_Description 

Confining 

Strain 

Rate.  6 

Temperature 

U  «  Unfrozen 

CIUC 

0.1,  2,  5 
and  10 

3  X  10-» 
(moderate) 

Rx)om 

FI  5=  Frozen 

UUC 

As  above 

As  above 

-10*  C 

F2  =  Frozen 

UUC 

0.1  &  5 

10-«(slow) 

lO'S(fast) 

-.5*C 
-20*  C 

FC  =  Frozen 

CIUC 

2&  10 

3  X  10‘5 

-10*  C 

B.  TEST  SYSTEMS 

1)  Manchester  Fine  Sand  (MFS)  at  55  &  95%  relative  density,  plus  three  ice-sand 
mixtures. 

2)  Boston  Blue  Clay  at  overconsolidation  ratio  =  4.0. 

3)  Fairbanks  Silt  at  natural  water  content. 

4)  Glass  Beads  at  95%  relative  density,  coarse  and  fine  silt  size,  plus  one  non— wetting. 


12 


PEAK 


OTTAWA  SANO  20-30 
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T  -  -7.fi*C 
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SOIL  STRENGTH 
♦  OILATANCY~ 

CXPECTEO  SLOPE 

CHANCE  V 


SOIL 

STRENGTH 
^  STRUCTURAL 
HINDRANCE' 


ICE  STRENGTHENING 


ICE  STRENGTH 


HYPOTHETICAL 

DIVISION 


VOLUME  FRACTION  OF  SAND  V,/V, 


(b)  Proposed  Mechanism  Map  for  Unconfined  Compressive  Strength  of  Frozen 
Ottawa  Sand  20—30  (Gonghnonr  and  Andersland,  1968)  at  —7.6*  C  and  Applied 

Strain  Bate  of  4.4x10“^  s“^ 


Figure  1.1: 


Proposed  Structure  and  Strength  Mechanisms  for 
Frozen  Sand  from  Ting  et  al.  (1983) 
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2.  EXPERIMENTAL  PROCEDURES 


2.1  INTRODUCTION 

A  ngnificant  effort  (approximately  one  year)  went  into  development  of  high 
pressure  triaxial  testing  equipment  for  both  the  frozen  and  unfrozen  testing  programs. 
Principal  attributes  of  the  frozen  testing  include:  lubricated  end  platens  (via  ice  caps); 
on-specimen  axial  strain  measurements  (precision  of  about  0.002%);  volumetric  strains 
(precision  of  ±0.2%)  from  volume  of  cell  fluid;  internal  load  cdl  to  measure  axial  force;  and 
good  temperature  control  (top  of  specimen  at  -9.4*  C,  bottom  at  -9.7*  C  and  fluctuations  of 
±0.05*  C  for  majority  of  frozen  tests).  In  addition,  high  pressure  triaxial  testing  capabilities 
(including  lubricated  end  platens  and  on— specimen  axial  strain  measurements)  were 
developed  for  characterizing  the  unfrozen  behavior  of  MFS  as  a  function  of  density  and 
consolidation  (confining)  pressure.  This  section  describes  the  test  material  and  the 
laboratory  facilities,  testing  equipment  and  experimental  procedures  used  in  testing  frozen 
and  unfrozen  Manchester  Fine  Sand. 

2.2  DESCRIPTION  OF  MANCHESTER  FINE  SAND 

The  sand  used  in  this  testing  program,  known  as  Manchester  Fine  Sand  (MFS),  was 
collected  at  the  Plourde  Sand  and  Gravel  Company  located  in  Hooksett,  NH  approximatdy 
10  miles  north  of  Manchester  along  the  banks  of  the  Merrimack  River.  This  sand  is  a  river 
bed  deposit  with  nearly  horizontal  bedding  planes.  The  gradation  can  vary  widely  between 
adjacent  bedding  planes  and  the  sand  for  this  testing  series  was  collected  from  several  of 
these  with  an  efibrt  to  obtain  representative  material. 

An  initial  mechanical  sieve  analysis  of  the  collected  material  showed  that  the  fines 
content,  i.e.,  the  material  finer  than  0.074mm,  was  approximately  20%.  The  sand  was 
subsequently  processed  by  dry  sieving  with  a  #200  sieve  to  remove  a  portion  of  the  fines. 
The  processed  batch  was  thoroughly  remixed  and  subsequent  sieve  analyses  indicated  that 
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the  fines  content  averaged  approximately  7%  (see  Fig.  2.1).  The  coefficient  of  uniformity 
for  this  sand  is  approximately  2.4  and  the  coefficient  of  concavity  is  approximately  1.1, 
based  on  the  average  of  four  sieves.  The  material  passing  the  No.  200  sieve  (fines)  is 
nonplastic.  The  Manchester  Fine  Sand  tested  in  this  research  is  classified  as  SP-SM,  a 
poorly  graded  fine  sand  with  7%  nonplastic  silt  size  particles  according  to  the  Unified  Soil 
Classification  System. 

Other  pertinent  characteristics  of  Manchester  Fine  Sand  are: 

•  Mainly  quartz  and  feldspar  (with  some  mica  flakes)  having  sub-angular 
particles  (Martin  et  al.  1981). 

•  Specific  gravity  =  2.688  ^  0.003. 

•  Minimum  dry  density  =  1408  kgf/m^;  maximum  dry  density  =  1701  kgf/m^. 

•  See  Fig.  2.2  for  unfrozen  water  content  versus  temperature. 

•  Although  MFS  is  significantly  more  compressible  than  Standard  Ottawa  Sand, 
it  shows  little  evidence  of  particle  crushing  at  a' c  =  10  MPA  (from  grain  size 
distribution  and  SEM  analyses). 

2.3  COLD  ROOM  FACILITY 

The  Low— Temperature  Testing  Fadlity  was  completed  in  1988  with  a  combination 
of  grants  from  NSF  and  two  MIT  projects,  ARO/PACT  and  the  Center  for  Scientific 
Excellence  in  Offshore  Engineering.  Constructed  and  located  inside  a  regular 
room— temperature  laboratory,  the  facility  has  three  separate  rooms  with  independent 
temperature  control:  the  v^tibule;  the  growth  room;  and  the  testing  room  (Fig.  2.3).  The 
vestibule  is  used  as  a  temperature  buffer  between  the  outside  laboratory  and  the  inner  cold 
rooms.  It  cuts  down  on  the  amount  of  humidity  that  reaches  the  two  inner  rooms  and  also 
serves  as  a  staging  area  for  preparing  and  trimming  frozen  soil  specimens.  It  normally  has 
a  temperature  of  —4  to  — 5*C.  The  growth  room,  located  to  the  rear  of  the  vestibule  and 
maintained  at  a  temperature  of  0*C,  is  used  for  specimen  freezing.  The  testing  room 
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contains  the  triaxial  apparatus,  loading  frame  and  the  environmental  chamber.  It  can  be 
regulated  to  —40*  C,  and  is  normally  maintained  at  a  temperature  a  few  degrees  colder  than 
the  desired  set  point  temperature  for  a  particular  test.  Temperature  control  of  the  triaxial 
cell  is  accomplished  by  heating  a  small  environmental  chamber  (Fig.  2.3)  located  inside  the 
testing  room. 

Each  room  is  composed  of  sheet  metal  walls  with  styrofoam  insulation.  The  testing 
room  has  a  double-paned  glass  window  to  allow  for  observation  of  the  testing  equipment 
from  the  outside,  room— temperature  laboratory.  This  window  is  heated  to  prevent  idng. 
The  testing  and  growth  rooms  are  fitted  with  feed-through  ports  to  allow  for  the  passage 
of  instrumentation  cables  and  power  lines.  These  are  used  to  connect  a  thermoregulator 
system,  load  cell,  displacement  transducers,  pressure  transducer  and  temperature  sensors  to 
the  controller  circuits,  power  supplies  and  data  acquisition  system  located  outside  the 
harsh  environment  of  the  testing  room.  See  Andersen  (1991)  for  a  comprehensive 
description  of  these  systems  and  devices.  Only  the  actual  sensing  devices,  the  circulating 
fan  and  heat  source  of  the  thermoregulator  system  and  the  pressure/volume  controller 
system  are  located  inside  the  cold  room. 

The  temperature  inside  the  testing  room  fluctuates  about  a  control  point.  The  air 
temperature  fluctuates  by  ±1.5’ C  about  the  desired  set  point  and  has  a  typical  period  of 
400  to  600  seconds.  The  testing  room  has  programmable  defrost  cycles  which  are  required 
for  continuous  operation.  These  defrost  cycles  produce  temperature  spikes  of  about  3’C 
above  the  normal  cyclic  oscillation  and  can  last  as  long  as  20  minutes. 

Testing  is  performed  inside  an  environmental  chamber  located  within  the  testing 
room.  The  environmental  chamber  is  a  double— walled  container  composed  of  sheets  of 
styrofoam  mounted  on  a  wooden  frame.  Figure  2.4  shows  a  schematic  of  this  chamber  and 
the  location  of  all  the  devices  within  it.  The  front  door  of  the  chamber  is  constructed  of 
two  sheets  of  plexiglas  mounted  on  an  aluminum  frame  to  allow  for  direct  observation  of 
the  testing  equipment  during  operation.  The  chamber  is  mounted  on  top  of  the  loading 
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frame  and  completely  encloses  the  triaxial  cell,  the  loading  platen,  the  crosshead  of  the 
loading  frame  and  the  oil  reservoir  for  the  pressure/volume  controller. 

The  temperature .  inside  of  the  environmental  chamber  is  controlled  by  a 
thermoregulator  system  composed  of  a  continuously  circulating  fan  (flow  rate  =  45cfrn),  a 
heat  source  mounted  inside  of  a  metal  box  and  a  mercury  contact  switch  preset  at  the 
desired  testing  temperature.  The  system  operates  by  adding  increments  of  heat  to  the 
chamber  as  heat  is  lost  to  the  colder  testing  room.  The  temperature  in  the  chamber 
fluctuates  around  the  desired  set  point  by  ±0.15*0  and  has  a  slight  gradient  of  about 
0.2*  C  (warmer  near  the  base). 

A  new  environmental  chamber  was  constructed  during  the  summer  of  1991; 
however,  most  of  the  aforementioned  characteristics  were  retained.  The  main  purpose  for 
the  new  chamber  was  to  create  a  more  stable  temperature  of  the  air  surrounding  the 
triaxial  cell.  For  example,  a  second  heat  source  was  added  to  improve  temperature 
stability.  A  similiar  environmental  chamber  was  also  constructed  for  the  unfrozen  soil 
testing  that  occurs  in  the  regular  room-temperature  laboratory. 

Digital  data  acquisition  is  an  integral  part  of  the  Low  Temperature  Testing  Facility. 
Output  voltages  from  transducers  and  thermistors  from  the  frozen  and  unfrozen  testing 
apparati  are  coimected  to  the  Geotechnical  Laboratory’s  Central  Control  Acquisition 
System.  The  system  has  a  140  channel  capability  and  connects  to  every  room  in  the 
laboratory  complex.  The  system  is  designed  around  the  Hewlett  Packard  3477A  data 
acquisition/control  unit.  This  is  a  very  low  noise  integrating  analog  to  digital  converter. 
It  has  17  bit  precision  and  auto  ranging  capability  resulting  in  one  micro  volt  sensitivity 
and  a  100  volt  range.  Any  type  of  direct  current  transducer  can  be  connected  to  this 
network  without  the  uncertainty  of  additional  signal  conditioners.  The  Low  Temperature 
Testing  Facility  currently  uses  25  data  acquisition  channds. 
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2.4  TRIAXIAL  TESTING  EQUIPMENT 

The  triaxial  cells  used  for  the  testing  program  were  manufactured  by  Wykeham 
Farrance  Engineering  Ltd.  of  Slough,  England.  They  were  first  used  at  MIT  25  years  ago 
for  investigating  the  strength-deformation  properties  of  stabilized  soils  (sponsored  by  the 
Army  Material  Command).  The  cells  were  modified  to  accommodate  lubricated  end 
platens,  internal  force  measurement,  on-spedmen  axial  strain  measurement  and  internal 
temperature  measurement. 

Figure  2.5  is  a  side  view  of  the  triaxial  cdl.  As  part  of  the  cells’  modification,  the 
bottom  pedestal  and  top  cap  were  replaced  by  enlarged  "lubricated"  platens  that  could 
accommodate  radial  deformation  of  the  soil  spedmen  during  shear.  The  spedmens  initially 
measure  1.4  inch  diameter  and  the  end  platens  are  1.6  inch  diameter.  The  modified  top 
cap  is  a  floating  top  cap,  meaning  that  it  is  not  rigidly  connected  to  the  loading  piston. 
The  cylinder  wall  of  the  triaxial  cdl  has  been  replaced  with  a  much  longer  one  to 
accommodate  the  length  of  the  new  top  cap  and  internal  load  cell.  The  modified  triaxial 
cell  used  for  unfrozen  tests  is  essentially  identical,  except  drainage  connections  have  been 
provided  to  both  the  bottom  and  top  of  the  spedmen. 

Silicone  510  oil  manufactured  by  the  Dow-Corning  Company  is  used  in  the  triaxial 
cell  as  the  pressurizing  medium.  It  was  sdected  for  its  low  viscosity  at  low  temperatures, 
its  nonconducting  properties  (necessary  for  internal  force  and  temperature  measurements 
with  dectronic  devices)  and  its  compatibility  with  the  prophylactic  membranes  used  to 
surround  the  frozen  and  unfrozen  sand  spedmens.  The  manufacturer  recommended  its  use 
for  temperatures  down  to  -51*  C.  Temperature  measurements  inside  the  triaxial  cell 
during  frozen  tests  indicate  that  the  silicone  oil  acts  as  an  effident  temperature  buffer. 
Data  taken  during  individual  tests  show  a  temperature  variation  of  ±0.05*0  at  any 
individual  location.  While  there  is  no  obvious  cyclic  component  to  the  variation,  there  is  a 
0.3*  C  temperature  gradient  within  the  cell.  This  gradient  is  opposite  (warmer  near  the 
top)  to  the  temperature  gradient  measured  in  the  environmental  chamber.  It  is  most  likdy 
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caused  by  heat  generated  by  the  internal  instrumentation.  We  are  currently  exploring 
methods  for  internal  circulation  of  the  silicone  oil  to  eliminate  the  temperature  gradient. 

The  cell  pressure  (and  the  back  pressure  in  the  unfrozen  tests)  is  applied  using  a  DC— 
servomotor  based  pressure/volume  control  system  which  is  capable  of  continuously  sensing 
and  adjusting  the  fluid  pressure  to  maintain  a  target  value.  A  DC-servomotor  is  attached 
to  a  ball-screw  actuator  which  drives  a  piston  into  a  reservoir  of  silicon  oil  (Fig.  2.6).  The 
motion  of  the  piston  then  serves  to  increase  or  decrease  the  fluid  pressure.  A  closed  loop 
analog  circuit  is  used  to  control  the  motor  for  the  frozen  tests.  This  analog  system  was  also 
initially  used  to  control  the  cell  pressure  in  the  unfrozen  tests  but  now  a  digital  closed  loop 
circuit  is  employed.  This  newer  circuit  system  provides  for  complete  computer  automation 
of  three  axes  (axial,  cell  and  pore  pressures)  during  the  testing  process  (see  next 
paragraph). 

The  current  frozen  soil  testing  system  uses  a  Wykeham  Farrance  T-n57  screw  driven 
loading  frame  having  a  nominal  capacity  of  10,000  pounds.  An  AC  motor  and  gear  box 
advance  the  loading  platen  at  a  constant  rate  of  deformation.  Due  to  the  flexibility  of  the 
loading  frame  and  the  use  of  lubricated  end  platens,  the  actual  rate  of  deformation 
experienced  by  the  specimen  is  not  constant,  but  increases  slowly  and  reaches  its  maximum 
at  about  0.5%  strain.  The  initial  strain  rate  is  roughly  one  order  of  magnitude  less  than 
the  maximum  or  nominal  strain  rate.  Thus  the  tests  are  not  really  sheared  at  a  constant 
strain  rate  until  reaching  the  "large  strain"  region.  (Note:  a  new  two  axis  computer  control 
system  vdll  be  implemented  to  decrease  the  strain  rate  variability  in  future  tests).  The 
unfrozen  soil  test  system  initially  used  a  Wykeham  Farrance  screw  driven  loading  frame. 
This  mechanical  system  was  replaced  by  a  hydraulic  loading  jack  capable  of  loads 
exceeding  8000  pounds.  The  new  hydraulic  system,  along  with  the  aforementioned 
DC-servomotor  based,  pressure  control  system,  allows  for  complete  computer  automation 
of  the  unfrozen  testing  program. 
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The  deviator  force  applied  to  the  specimen  is  measured  by  a  load  cell  located  inside 
the  triaxial  chamber  at  the  end  of  the  loading  piston.  This  avoids  measurement  errors  due 
to  friction  between  the  loading  piston  and  its  brass  bushing  guide.  The  load  cell  contacts 
the  floating  top  cap  with  a  cylindrical  seating  piece.  This  piece  is  in  contact  with  the 
entire  flat  ^d  of  the  cylindrical  top  cap  and  has  a  small  circular  lip  with  an  inside  taper  to 
provide  a  guide  during  set  up  when  the  load  cdl  and  piston  must  be  lowered  onto  the  top 
cap.  This  lip  also  restricts  lateral  motion  of  the  top  cap  during  shear. 

Axial  strains  are  measured  by  two  Alternating  Current  Differential  Transformers 
(ACDT’s)  mounted  in  diametrically  opposite  positions  around  the  circumference  of  the 
sand  specimen  on  a  specially  designed  yoke  (Fig.  2.7).  The  upper  piece  holds  the  barrels  of 
the  displacement  transducers  and  the  cores  rest  on  the  lower  yoke.  Each  yoke  is  spring- 
loaded  and  contacts  the  specimen  at  third  points  around  the  circumference.  The  face  of 
each  contact  foot  is  a  flat  rectangle  approximately  0.125  inches  wide  and  0.25  inches  long. 
During  set  up,  three  alignment  rods  run  between  the  upper  and  lower  yoke.  These  are  held 
in  position  with  0.063  inch  diameter  hardened  steel  pins.  The  function  of  these  rods  is  to 
aid  in  the  placement  of  the  yoke  on  the  specimen  and  to  fix  the  appropriate  gage  length. 
After  placing  the  yoke  on  the  specimen,  these  small  pins  are  removed  and  the  alignment 
rods  drop  freely  away  from  the  yoke  allowing  it  to  be  free  floating  and  move  unimpeded  as 
the  specimen  deforms. 

Interchangeable  cell  pressure  (and  pore  pressure  for  the  unfrozen  tests)  transducers 
are  mounted  directly  on  the  triaxial  cell.  These  have  capacities  of  200  psi,  2000  psi  or  5000 
psi  depending  on  the  particular  test  requirements. 

Volume  changes  of  the  frozen  sand  specimens  during  shear  were  computed  from  the 
motions  of  the  cell  pressure/volume  controller  piston  and  the  triaxial  cell  loading  piston. 
Displacements  of  these  are  measured  by  DC— DC  type  Linear  Variable  Differential 
Transformers  (LVDT’s)  mounted  on  each  piston.  The  corresponding  volume  changes  are 
computed  from  the  cross-sectional  area  of  each  piston.  Various  "corrections”  were  then 
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applied  to  these  meastuements  in  order  to  account  for  leakage  of  the  cell  fluid, 
compressibility  of  the  cell  fluid  and  triaxial  chamber  and  flexure  of  the  base  of  the  triaxial 
cell  (Chapter  3  of  Andersen  1991).  This  method  of  measurement  results  in  a  maximum 
volumetric  strain  error  of  *0.20%.  The  error  increases  with  the  test  duration  (i.e.,  the  slow 
tests  are  subject  to  the  largest  error). 

Volume  changes  during  drained  shear  tests  on  unfrozen  sand  are  measured  from 
motions  of  the  back  pressure  controller  piston.  This  is  a  direct  measure  of  the  fluid  passing 
into  or  out  of  the  specimen  and  is  not  subject  to  the  same  complications  as  for  the  frozen 
tests. 

2.5  FROZEN  MANCHESTER  FINE  SAND  TEST  PROCEDURES 

2.5.1  Specimen  Preparation 

The  molds  used  for  the  preparation  and  freezing  of  the  Manchester  Fine  Sand 
specimens  were  obtained  from  CRUEL  and  were  the  same  molds  used  previously  by  Martin 
et  al.  (1981).  Figure  2.8  presents  a  series  of  schematic  drawings  of  one  of  these  molds  in 
different  stages  of  operation.  Figure  2.8(a)  is  the  mold  prepared  for  the  wet  tamping 
procedure,  Fig.  2.8(b)  is  the  mold  prepared  for  multiple  sieve  pluviation.  Fig.  2.8(c)  is  the 
mold  prepared  for  saturation  and  Fig.  2.8(d)  is  the  mold  prepared  for  freezing. 

The  vast  majority  of  the  spedm^s  were  prepared  using  multiple  sieve 
pluviation  (MSP)  as  outlined  by  Miura  and  Told  (1982).  This  produced  spedmens  of 
uniform  density  with  relative  densities  ranging  from  20  to  100%.  In  this  procedure  (set-up 
shown  in  Fig.  2.9),  oven  dried  sand  is  poured  into  a  funnel  mounted  in  a  frame  above  the 
mold.  The  sand  falls  out  of  the  funnel,  through  a  vertical  plexiglas  tube  and  across  four  3 
inch  diameter  sieves  (two  with  No.  20  mesh  sizes  and  two  with  No.  10  mesh  sizes).  The 
sieves  act  to  disperse  the  raining  doud  of  sand.  The  sand  then  falls  across  a  knife  edge 
guide  and  into  the  mold. 
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After  pluviation,  the  dry  sand  spedm^  are  taken  into  the  growth  room 
maintained  at  0*  C.  The  upper  side  of  each  mold  is  connected  to  one  reservdr  and  the 
lower  side  is  connected  to  another  reservoir.  A  28-^nch  Hg  (95  kPa)  vacuum  is  drawn  on 
the  specimens  and  deaired,  deionized,  distilled  water  is  passed  through  the  molds  from  the 
base  to  the  top  by  gravity  feed  under  an  initial  head  of  approximately  five  feet  of  water. 
After  several  pore  volumes  of  water  pass  through  the  specimens,  the  upper  and  lower 
reservoirs  are  raised  to  the  same  elevation  above  the  specimens.  This  is  maintained  under 
vacuum  over  night. 

The  following  morning,  the  vacuum  is  slowly  removed,  the  upper  reservoir  is 
disconnected  and  the  lower  reservoir  is  placed  at  the  elevation  of  the  top  of  the  specimens. 
The  brass  top  caps  of  the  specimen  preparation  molds  are  removed  and  replaced  with 
freezing  caps.  The  filter  papers  are  retained  on  the  top  of  each  specimen.  The  freezing 
caps  are  placed  in  direct  contact  with  the  filter  paper.  Refer  to  Figure  2.8(d).  Antifreeze 
fluid  at  -15*  C  is  circulated  through  the  freezing  caps  and  the  specimens  are  frozen  from 
the  top  down,  while  allowing  for  drainage  through  the  base  of  the  molds  and  into  the  lower 
reservoir.  The  entire  freezing  process  lasts  approximately  four  hours.  After  freezing,  the 
specimen  (with  surrounding  plexiglas  split  sleeve)  is  removed  from  the  mold  and  stored  in  a 
freezer  until  trimmed  for  testing. 

The  actual  trimming  of  the  specimen  takes  place  in  the  vestibule  room  of  the 
cold  room  facility.  This  room  is  maintained  at  -4  to  -5*  C.  The  specimens  are  taken  from 
the  testing  room  at  —15*  C  to  the  vestibule  for  a  maximum  of  5  minutes  for  trimming.  This 
prevents  melting  due  to  handling.  One  end  of  the  specimen  is  extruded  approximately  0.25 
inches  out  of  the  split  sleeve  and  a  hose  clamp  is  placed  on  the  split  sleeve  to  hold  the 
specimen  firmly  in  position.  The  trimming  is  done  by  hand  with  sharpened  sted  knives. 
After  trimming  one  end,  the  specimen  is  partially  extruded  from  the  other  end  and  the 
process  is  repeated.  In  order  to  measure  the  parallelism  of  the  ends,  the  specimen  is  taken 
out  of  the  split  sleeve  and  placed  on  a  metal  stand  with  a  spring  loaded  depth  gage 


sensitive  to  0.0001  inch.  The  height  is  selectively  trimmed  to  within  a  0.001  inch  tolerance 
across  the  diameter  of  the  specimen. 

The  diameter  of  the  specimen  is  measured  with  a  hand  held  micrometer  with  a 
direct  readout  to  0.001  in.  at  eight  locations  along  the  length  (at  each  end  and  third  points, 
two  measurements  at  90*  to  each  other  at  each  point).  The  height  of  the  spedmen  is 
measured  at  four  locations  spaced  at  roughly  90*  around  the  circumference.  The  averages 
of  the  diameter  and  height  measures  are  used  subsequently  to  compute  the  specimen 
volume.  The  entire  trimmed  specimen  is  then  weighed  on  an  electronic  scale  sensitive  to 

±0.01  g. 

2.5.2  Specimen  Set-Up 

For  those  specimens  tested  with  "lubricated”  ends,  an  ice  cap  was  frozen  onto 
each  end  and  subsequently  trimmed  in  the  same  manner  as  the  frozen  sand.  In  order  to 
properly  align  the  frozen  sand  specimen  between  the  base  pedestal  and  floating  top  cap  of 
the  triaxial  cell,  0.063  in.  diameter  holes  are  drilled  through  the  center  of  each  ice  cap  and 
into  the  frozen  sand  using  a  drill  bit  mounted  in  a  plexiglas  boring  tool  which  fits  around 
the  outside  of  the  plexiglas  split  sleeve.  These  holes  are  bored  approximately  0.125  in.  deep 
and  mate  with  0.063  in.  diameter  hardened  sted  pins  mounted  in  the  center  of  the  base 
pedestal  and  top  cap. 

Prior  to  placement  of  the  trimmed  frozen  sand  specimen  in  the  triaxial  cell,  the 
base  pedestal  and  top  cap  are  covered  with  high— vacuum  silicone  grease  and  a  prophylactic 
membrane  is  placed  on  the  base  pedestal  with  two  0-^ngs.  The  specimen  is  placed  on  the 
base  pedestal  (bottom  end  down)  and  the  top  cap  is  lowered  onto  the  upper  surface  of  the 
specimen  using  a  specially  designed  alignment  jig  mounted  directly  on  the  base  of  the 
triaxial  cell.  An  axial  force  is  applied  by  hand  to  the  specimen  through  the  alignment  jig 
in  order  to  force  the  excess  grease  from  between  the  specimen  and  the  brass  surfaces  of  the 
base  pedestal  and  top  cap  and  to  bring  the  frozen  sand  specimen  into  vertical  alignment. 
The  grease  also  acts  to  hold  the  specimen  and  top  cap  on  the  base  pedestal  during  the  set 
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up  procedure.  The  extruded  grease  is  wiped  off  and  the  rubber  membrane  is  stretched  from 
the  base  pedestal  over  the  specimen  and  onto  the  top  cap.  One  of  the  0-rings  on  the  base 
pedestal  is  raised  to  the  top  cap  to  secure  the  membrane  in  place. 

The  axial  displacement  yoke  is  placed  on  the  outside  of  the  membrane  at 
roughly  the  third  points  along  the  specimen  so  that  the  active  gage  length  includes  the 
central  portion  of  the  specimen.  After  placement  of  the  yoke  on  the  specimen  and  removal 
of  the  alignment  posts,  the  gage  length  is  measured  using  an  optical  scope  with  a  veneer 
that  has  a  direct  read  out  of  0.001  of  an  inch. 

With  the  axial  displacement  yoke  in  position,  the  top  of  the  triaxial  cell  is 
lowered  onto  the  base  using  a  guide  stand  to  prevent  accidental  disturbance  of  the  axial 
displacement  yokes.  The  load  cell  is  gently  lowered  onto  the  top  cap  and  the  top  and  base 
of  the  triaxial  cell  are  bolted  together.  The  loading  piston  is  then  locked  into  position  and 
the  entire  triaxial  cell  is  placed  in  the  loading  frame.  Each  of  the  electronic  devices  are 
connected  to  the  power  supplies  and  data  acquisition  system  and  the  output  level  of  each 
device  is  checked  to  verify  that  it  was  within  the  normal  range. 

A  slight  preload  (0.1  to  0.3  MPa)  is  applied  to  the  specimen  through  a  manual 
crank  on  the  loading  frame.  It  is  reapplied  at  various  times  throughout  the  preshear  stage 
of  the  test  because  it  gradually  reduces  over  time  due  to  squeezing  out  of  the  grease  and 
creep  of  the  ice  cap  and  frozen  sand. 

2.5.3  Specimen  Shearing 

The  triaxial  cell  is  filled  with  silicone  oil  using  a  diaphragm  pressure  reservoir 
located  inside  the  cold  room.  The  cell  pressure  is  applied  and  the  specimen  is  allowed  to  sit 
in  the  triaxial  cell  under  pressure  with  the  preload  until  both  the  temperature  and  the 
leakage  rate  of  oil  out  of  the  triaxial  cdl  stabilize.  This  presheai  pressurization  phase  lasts 
longer  than  the  expected  shearing  time  of  each  test. 

After  the  preshear  period,  the  compressibility  of  the  cell  fluid  is  measured  by 
cycling  the  cell  pressure  and  measuring  the  corresponding  volume  change.  Just  prior  to 
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shearing,  the  small  preload  is  removed  and  the  zero  of  the  load  cdl  is  recorded.  The 
deviator  stress  is  applied  to  the  spedmen  with  the  screw  driven  loading  frame.  In  general, 
the  specimens  are  sheared  to  axial  strains  in  excess  of  20%.  After  shearing,  the  cell 
pressure  is  relieved  and  the  triaxial  cell  is  removed  from  the  loading  frame  and 
disassembled.  The  specimen  is  photographed  and  then  removed  from  the  base  pedestal. 
The  final  diameter  is  measured  near  the  ends  and  at  third  points  along  the  length. 

Calculation  of  the  dry  density  of  the  test  specimen  is  based  on  the  volume  of  the 
trimmed  specimen  before  shear  and  dry  weights  measured  after  pluviation  (minus  the 
weight  of  trimmings),  after  shearing  (which  may  contain  small  amounts  of  silicone  grease 
or  oil),  and  after  washing  to  remove  any  silicone  (which  may  also  lose  some  soil  particles). 
The  estimated  standard  deviation  in  the  reported  values  of  relative  density  is  ^2%.  Table 
A.l  presents  the  physical  properties  of  the  specimens. 

2.5.4  Precision  of  Stress— Strain  Parameters 

Chapter  4  of  Andersen  (1991)  presents  an  "error  analysis"  of  results  from  nine 
pairs  of  tests  in  order  to  estimate  the  likely  precision  of  the  reported  test  data.  The 
following  summarizes  the  principal  conclusions,  where  COY  is  the  coefficient  of  variation 
and  SD  is  the  standard  deviation. 

1)  For  Young’s  modulus  wherein  the  two  ACDT’s  showed  "good"  or  "fair"  (as 
opposed  to  "poor")  agreement: 

COY  =  9.5%  with  SD  =  ±4.5%. 

2)  For  the  upper  yield  stress  and  the  conesponding  axial  strain: 

COY  =  7%  with  SD  =  ±3.5%  for  stress 

COY  =  8%  with  SD  =  ±7%  for  strain. 

3)  For  the  peak  strength  and  the  corresponding  axial  strain  from  tests  having 
"good",  "fair"  and  "poor"  (as  opposed  to  "very  poor")  stability.  These 
qualifiers  regard  the  lateral  sliding  stability  at  large  strains  due  to  the 
lubricated  end  conditions: 
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COV  s  3.5%  with  SD  =  *2.5%  for  stress 
COV  =  13%  ¥rith  SD  s  ^^9%  for  strain. 

Andersen  (1991)  also  describes  the  procedure  used  to  compute  the  instantaneous  cross- 
sectional  area  of  the  specimens  during  shear.  This  made  it  possible  to  correct  for  the 
effects  of  bulging  within  the  central  portion  or  flaring  at  the  ends  that  occurred  in  many  of 
the  tests. 

2.6  UNFROZEN  MANCHESTER  FINE  SAND  TEST  PROCEDURES 

Many  of  the  same  procedures  used  in  setting-up  frozen  specimens  are  also  used  on 
unfrozen  specimens.  The  major  difference  from  frozen  test  specimens  is  that  unfrozen 
specimens  are  formed  directly  on  the  pedestal  of  the  triaxial  cell  base.  It  is  important  to 
set-up  the  specimen  on  a  dean,  dry  system.  The  base  pedestal  and  top  cap  are  prepared 
by  first  flushing  all  drainage  lines  and  ports  with  water  to  remove  any  debris  (sand,  silicon 
oil,  etc.)  in  the  lines.  The  lines  are  then  flushed  with  pressurized  air  to  remove  the  water 
from  the  lines  and  dry  the  drainage  system.  Porous  stones  are  placed  in  the  drainage  ports 
in  the  pedestal  and  top  cap.  Lubricated  ends,  ranging  from  greased  latex  rubber  to  a 
silicone  coating,  are  then  applied  to  the  pedestal  and  top  cap.  A  thick  (0.014-inch)  rubber 
membrane  is  attached  to  the  pedestal  using  0— rings  and  vacuum  grease.  A  membrane 
expander  is  then  placed  over  the  pedestal,  and  a  slight  vacuum  (8-inch  Hg)  is  placed  on 
the  membrane  to  create  a  cylindrical  (1.4  inch  diameter)  volume  that  forms  the  specimen. 
Multiple  sieve  pluviation  is  used  to  form  specimens  (see  Fig.  2.10  and  Section  2.5.1),  after 
which  excess  sand  is  removed  from  the  top  of  the  newly— formed  specimen  and  the  top  cap 
is  placed  on  the  specimen  using  the  alignment  pg. 

An  internal  drainage  line  is  attached  to  the  top  cap  to  create  a  closed  system.  The 
small  vacuum  being  maintained  on  the  outside  of  the  membrane  (to  keep  it  open  during 
sand  placement)  is  now  transferred  to  the  specimen  and  increased  to  15-inch  Hg.  This 
vacuum  will  be  maintained  until  the  specimen  is  placed  under  an  equivalent  cell  pressure. 
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The  membrane  expander/sample  former  is  removed,  and  sample  dimensions  (height  and 
diameter)  are  measured  using  the  optical  scope.  Two  ACDTs  are  placed  on  the  specimen 
using  an  axial  displacement  yoke  assembly  similar  to  that  used  in  frozen  MFS  tests.  The 
gage  length  is  then  measured  using  the  optical  scope. 

The  top  of  the  triaxial  cell  is  lowered  onto  the  base  and  load  cell  is  gently  lowered 
onto  the  specimen.  The  top  and  base  of  the  cell  are  then  bolted  together  and  the  piston  is 
locked  in  place.  The  entire  cell  is  moved  over  to  the  loading  frame,  and  electronic  devices 
are  connected  and  checked.  The  cell  is  raised  so  that  the  piston  is  in  contact  with  the 
loading  point  of  the  load  frame.  The  cell  is  then  filled  with  oil. 

The  following  testing  procedures  are  those  used  with  the  computer  automated 
system.  Note  that  prior  to  computer  automation,  the  raising  and  lowering  of  cell  and  back 
pressures  was  performed  using  an  analog  circuit  loop  and  manual  control  techniques. 
Computer  automation  has  allowed  for  better  control  of  the  cell  and  back  pressure  during 
the  test.  A  cell  pressure,  equivalent  to  the  applied  vacuum  on  the  specimen,  is  placed  on 
the  specimen  to  provide  support.  As  the  cell  pressure  increases,  the  vacuum  is  slowly 
released  so  that  the  effective  stress  on  the  sample  remains  approximately  the  same.  When 
the  cell  pressure  is  reached  and  the  vacuum  has  been  totally  removed  from  the  spedmen, 
the  spedmen  is  flushed  with  carbon  dioxide  for  approximately  15  minutes  to  remove  as 
much  air  as  possible  from  the  pore  voids.  At  least  three  pore  volumes  of  deaired, 
de-4onized  water  are  then  passed  through  the  spedmen  (flowing  bottom  to  top)  to  displace 
the  carbon  dioxide. 

The  spedmen  is  saturated  by  increasing  the  cell  pressure  and  back  pressure,  in  equal 
increments  to  maintain  the  same  effective  stress,  until  a  cell  pressure  of  1.1  MPa  and 
backpressure  of  1.0  MPa  are  obtained.  The  spedmen  is  allowed  to  saturate  overnight 
under  these  conditions.  The  change  in  volume  of  the  spedmen  during  overnight  saturation 
is  then  inferred  from  overnight  volume  measurement  as  well  as  a  leakage  rate  from  the 
back  pressuring  system  under  constant  pressure  conditions. 
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The  following  day,  snccessfiil  satniation  is  confirmed  by  measuring  the  B— value 
(=An/A(rc)  using  computer  control.  With  a  successful  B— value  obtained  (i.e.,  at  least 
95%  within  2  minutes),  the  specimen  is  hydrostatically  consolidated  by  incrementally 
raising  the  cell  pressure  to  obtain  the  desired  preshear  effective  stress  (Note:  specimens 
sheared  at  an  effective  confining  pressure  of  0.1  MPa  do  not  undergo  this  consolidation 
process).  The  changes  in  pore  volume  are  used  to  compute  the  volumetric  strain  during 
consolidation. 

The  specimen  is  sheared  after  allowing  some  time  for  secondary  compression 
(usually  two  or  more  hours).  In  undrained  shear,  the  backpressure  drainage  lines  are  closed 
and  pore  pressures  developed  during  shear  are  measured.  In  drained  shear  tests,  the  back 
pressure  drainage  lines  remain  open  and  volume  changes  are  measured  as  the  controller’s 
piston  moves  in  or  out  to  maintain  the  constant  back  pressure.  The  specimen  is  usually 
sheared  to  axial  strains  in  excess  of  20  percent. 

After  shearing,  the  cell  and  back  pressures  are  relieved,  oil  is  removed  from  the  cell, 
and  the  cell  removed  from  the  loading  frame  and  disassembled.  The  specimen  is 
photographed  and  specimen  form  is  noted.  The  specimen  is  frozen  while  still  on  the 
pedestal,  then  the  frozen  specimen  is  removed  from  the  cell  base  for  measurements  of  soil 
weights  before  and  after  oven  drying.  These  data  are  used  to  compute  the  preshear  void 
ratio  and  relative  density  (after  accounting  for  volume  changes  for  drained  shear  tests). 
Similar  variations  in  specimen  properties  exist  for  the  unfrozen  MFS  specimens  as  noted  for 
frozen  MFS  specimens. 
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Figme  2.1:  Mechanical 
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Figoie  2.3:  Plan  View  of  MIT  Cold  Room  Facility 
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Figure  2.4: 


Schematic  View  of  Environmental  Chamber  and 
Triasdal  CeU  (Before  Addition  of  Second  Heat  Source) 
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Pressurizing  Piston 


Figure  2.6:  Schematic  Drawing  of  Mechanical  Portion  of 

Pressure/ Volume  Controller 
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Figure  2.8:  Schematic  Drawing  of  Specimen  Preparation  Molds 

(after  Martin  et  al.  1981) 
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3.  BEHAVIOR  OF  FROZEN  MANCHESTER  FINE  SAND 
IN  TRIAXIAL  COMPRESSION 


This  section  presents  the  results  of  the  triaxial  compression  tests  conducted  on 
frozen  Manchester  Fine  Sand.  Most  of  this  presentation  is  abstracted  from  Andersen’s 
(1991)  doctoral  thesis  that  focused  on  a  compr^ensive  evaluation  of  the  effects  of  relative 
density  (Dr),  stiwu  rate  (ir)  and  confining  pressure  (ac)  on  stress-strain  behavior  at  one 
temperature  (-9.5*  C).  After  outlining  the  scope  of  the  testing  program,  which  includes 
subsequent  tests  at  varying  temperature,  a  qualitative  overview  of  the  stress-strain 
behavior  and  definitions  of  the  parameters  used  to  describe  the  stress-strain  response 
follow.  A  quantitative  treatment  of  the  stress-strain  behavior  at  -9.5*  C  is  then  presented 
in  two  major  sections,  small  strain  behavior  and  large  strain  behavior.  The  small  strain 
behavior  is  considered  in  terms  of  the  Young’s  modulus,  the  proportional  limit  and  the 
yield  offset  stress  at  10  strain.  The  large  strain  behavior  is  considered  in  terms  of  the 
upper  yield  region,  the  strain  hardening  or  strain  softening  region,  the  peak  strength  region 
and  the  volumetric  behavior. 

Section  3  also  presents  results  from  tests  conducted  at  different  temperatures. 
These  data  and  their  evaluation  are  incomplete.  Section  3  concludes  with 
recommendations  for  further  research,  namely:  completion  of  testing  at  varying 
temperatures;  and  devdopment  of  new  experimental  techniques  in  order  to  conduct  tests 
on  sand  that  is  frozen  after  application  of  the  confining  pressure  rather  than  vice  versa. 
This  will  allow  for  the  evaluation  of  how  the  magnitude  of  the  pre-fireezing  effective  stress 
acting  on  the  sand  skeleton  affects  strength-deformation  properties. 

Appendix  A  contains  the  tabulated  data  used  to  construct  the  summary  plots 
presented  in  Section  3  (See  the  cover  page  of  Appendix  A  for  a  listing  of  these  tables). 
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. . . . '  -  WII.  i  ■■  III  .pjipjJ^JUIP 

I  3.1  SCOPE  OF  TESTING  PROGRAM 

Table  3.1  sommaiMes  the  scope  of  the  testing  program.  The  "principal"  program 
consists  of  tests  at  various  relative  densities,  confining  pressures  and  strain  rates.  The 
I  rdative  density  (Dr)  varied  from  20  to  100%  and  the  confining  pressure  (ctc)  varied  from 
0.1  to  10.0  MPa.  Three  nominal  strain  rates  (Ie)  were  used:  3xl0'>/sec,  3xlO’V8ec  and 
4xl0'</sec,  referred  to  as  "slow",  "moderate"  and  "fast"  tests,  respectively.  One  specimen 
.  was  tested  at  a  strain  rate  of  1.2xl0‘^/sec  (referred  to  as  "intermediate").  The  testing 
temperature  was  -9.55^0.3*  C,  except  for  one  test  having  a  temperature  of  about  -15*  C. 

I  Six  tests  are  presented  to  show  initial  results  from  the  new  temperature  testing 

program.  Temperatures  for  these  moderate  strain  rate  tests  range  from  —5  to  —28  *  C.  A 
I  confining  pressure  of  0.1  MPa  was  used  and  rdative  densities  range  from  41  to  48%. 

t 

3:2  DEFINITION  OF  PARAMETERS  AND  GENERAL  OVERVIEW  OF  STRESS- 
j  STRAIN  BEHAVIOR 

3.2.1  Definition  of  Parameters 

I  Figure  3.1  shows  the  graphical  construction  techniques  used  by  Andersen  (1991) 

I  to  obtain  various  stress— strain  parameters  from  the  triaxial  compression  test  program. 

The  figure  plots  engineering  axial  strain  (ca  =  AL/Lq)  versus  deviator  stress  (Q  =  ai  —  as) 

I  and  volumetric  strain  (cv  =  AV/Vq,  positive  for  dilation  =  expansion).  The  presentation 

I 

will  focus  on  the  following  parameters  and  behavioral  features: 
i  1)  Values  of  Young’s  modulus  (initial  slope  of  the  stress-strain  curve),  the 

proportional  limit  (point  on  the  stress-strain  curve  where  the  response  departs 
significantly  from  linear),  and  the  yield  of&et  stress  at  10'^  strain  (intersection  point 
of  the  stress-strain  curve  and  a  line  with  a  slope  of  the  Young’s  modulus  translated 
by  10*4  strain).  See  Fig.  3.1(a),  plotted  to  €a  ==  0.1%,  for  graphical  definitions. 

2)  Values  of  the  upper  yield  stress  (Quy),  the  knee  of  the  stress-strain  curve,  and  its 
corresponding  strain  (cy)  obtained  from  Fig.  3.1(b)  for  stress-strain  curves  plotted 
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to  6a  =  2.6%.  (Note:  although  plastic  strains  start  at  the  proportional  limit,  Fig. 
3.1(a),  they  only  become  really  significant  in  the  "upper  yield"  region). 

3)  The  nature  of  the  post  upper  yield  behavior,  i.e.,  the  degree  of  strain  hardening  or 
strain  softening  after  the  knee  as  per  Fig.  3.1(b).  Most  curves  exhibited  a  distinct 
double  yield  behavior  similar  to  that  of  fine-grained  polycrystalline  ice  (e.g.,  Mellor 
and  Cole  1982;  Cole  1987)  as  further  described  in  Section  3.2.2. 

4)  Values  of  the  peak  stress  (Qp)  and  corresponding  strain  (ep)  as  per  Fig.  3.1(c).  In 
general,  specimens  tested  at  low  confining  pressures  experienced  post  peak  strain 
softening,  wherein  those  at  the  highest  confining  pressures  did  not  (see  Section 
3.2.2). 

5)  The  volumetric  behavior,  including  the  rate  of  dilation  =  d£v/d6a,  as  per  Fig. 
3.1(d). 

3.2.2  Types  of  Stress-Strain  Curves 

The  general  shape  of  the  stress-strain  curves  as  a  function  of  the  relative 
density,  confining  pressure  and  strain  rate  can  be  classified  into  seven  different  types. 
Figure  3.2  presents  several  normalized  (i.e.,  deviator  stress  divided  by  the  upper  yield 
stress,  Q/Quy)  stress-strain  curves  that  serve  to  illustrate  visually  these  different  types  of 
behavior. 


A  qualitative  description  of  each  type  of  stress-strain  behavior  follows. 

Low  Confining  Pressure 

Type  I  The  peak  stren^h  occurs  at  the  upper  yield  point.  This  is  followed  by 
pronounced  strmn  softening. 

Type  n  The  upper  yield  stress  (knee  in  the  stress-strain  curve)  is  followed  by  a 
yield  ^op.  A  slight  amount  of  strain  hardening  occurs  after  this  yield 
drop  to  a  second  peak  having  a  strength  similar  to  the  upper  yield  stress. 
Some  strain  softening  occurs  at  very  large  strains. 

Type  m  The  upper  yield  stress  is  followed  by  a  plateau  in  the  stress-strain  curve. 

After  this  plateau,  there  is  a  significant  amount  of  strain  hardening  to  the 
peak  strength,  followed  by  strain  softening. 
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Type  IV  The  upper  yidd  ttxese  is  Mowed  by  immediate  strain  hardening  to  the 
peak  stren^.  After  the  peak,  there  is  a  considerable  amount  m  strain 
softening. 

For  the  low  confining  pressure  Type  I— IV  curves,  strain  softening  causes  the  large  strain 
post  peak  resistance  to  be  generally  less  than  the  upper  yield  stress. 

High  Confining  Pressure 

Type  V  The  upper  yidd  stress  is  followed  by  a  yidd  drop.  After  the  yidd  drop, 
there  is  a  slight  amount  of  strain  hardening.  The  deviator  stress  is  almost 
constant  to  very  large  strains. 

Type  VI  The  upper  yidd  stress  is  followed  by  a  small  plateau.  After  the  plateau, 
there  is  continuous  strain  hardening  to  a  peak  strength  at  very  large 
strains. 

Type  VU  The  upper  yidd  stress  is  followed  by  an  immediate  and  si^ficant  amount 
of  strain  hardening  to  the  peak  strength  at  very  large  strains. 

For  the  high  confining  pressure  Type  V— Vn  curves,  the  large  strain  resistance  will 
generally  be  larger  than  the  upper  yidd  stress. 

The  presentation  first  summarizes  Young’s  modulus,  etc.  to  represent  small 
strain  behavior,  and  then  various  aspects  of  large  strain  behavior  (upper  yidd  stress,  peak 
stress,  volumetric  strains,  etc.)  as  a  function  of  sand  density,  confining  pressure  and  strain 
rate. 


3.3  SMALL  STRAIN  BEHAVIOR 

The  small  strain  behavior  can  be  described  in  terms  of  the  measured  Young’s 
modulus,  the  proportional  limit,  and  the  yidd  ofiset  stress  at  10*^  strain. 

3.3.1  Youna’s  Modulus 

Figure  3.3  plots  Young’s  modulus  (E)  versus  rdative  density  (Dr)  for  all 
confining  pressures  and  strain  rates  from  tests  having  "good  to  fair"  agreement  between  the 
two  ACDTs  used  to  measure  the  on-spedmen  axial  strains.  Due  to  compliance  in  the 
triaxial  equipment,  the  initial  strain  rates  were  about  one— tenth  of  the  values  listed  in 
Table  3.1  (See  Table  A.2).  Evaluation  of  the  data  and  comparison  with  prior  results  for 
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ficozen  sand  show  the  following  trends  as  a  function  of  relative  density  (Dr),  confining 
pressure  (trc)  and  strain  rate  (c). 

Trend  in  Yonng*s  Modulus 


Variable 

This  Program 

Prior  Results 

Dr  =  25  to  95% 

No  clear 

Baker  and  Durfurst 

trend 

(1985):  Increase  of  20% 
mm  wave  velocity 
data  on  Ottawa  sand 

ffc  =  0.1  to 

15%  decrease 

No  rdiable  data 

10  MPa 

€  =  3  X  10‘*  to 

No  effect 

No  data 

4  X  lO’Vs 

Andersen  (1991)  used  the  isostrain  model  proposed  by  Counto  (1964),  which 
was  developed  to  predict  the  modulus  of  concrete,  in  order  to  model  the  small  strain 
behavior  of  frozen  MFS.  In  this  model,  the  composite  material  is  idealized  as  a  cube  of 
aggregate  (sand  grains)  embedded  in  a  cube  of  matrix  material  (ice). 

An  axial  force  is  applied  to  a  composite  cube  having  three  horizontal  sections 
(i.e.,  two  pure  ice  ends  and  a  central  section  having  an  isostrain  condition  imposed  between 
the  entire  volume  of  a  silicate  sand  prism  and  the  surrounding  ice)  and  the  corresponding 
deformation  of  each  section  is  computed  from  the  "known"  dastic  moduli  of  the  pore  ice 
and  the  sand  particles.  The  resulting  equation  for  the  composite  Young’s  modulus  (Ec) 
using  the  height  to  area  assumption  of  Counto  (1964)  is 


1 

E7 


1  -  = 

BTi 


1  -  rr. 


(3.1) 


Ei  +  Eg 


where  V,  is  the  volume  fraction  of  the  sand  particles,  Ei  is  the  Young’s  modulus  for  ice  and 
Eg  is  the  Young’s  modulus  of  the  sand  particles  (not  the  sand  skeleton).  If  the  geometry  of 
the  silicate  prism  in  the  central  section  is  that  of  a  perfect  cube,  a  similar  relationship 
predicts  values  of  Ec  about  20%  higher  than  Eq.  3.1. 
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Selecting  Es  =  90  GPa  and  Ei  =  7.5  GPa  for  sand  grains  and  ice,  respectively, 
Andersen  computed  Ec  =  25  ±  2.5  GPa  at  Dr  =  35%  and  Ec  =  29  3  GPa  at  Dr  =  100%. 

These  predictions  agree  very  weU  with  the  data  for  MFS  which  give  E  =  26.6  GPa  ±  4.5  SD 
and  hence  illustrate  the  value  of  considering  a  frozen  sand  system  as  a  composite  material. 
It  should  be  emphasized  that  the  much  higher  modulus  of  frozen  sand  compared  to  pure  ice 
does  not  result  from  the  stiffuess  of  the  sand  sJreieton,  which  is  much  less  than  that  of  ice, 
but  rather  from  the  very  high  modulus  of  sand  gridns  embedded  in  an  ice  matrix. 

3.3.3  Proportional  Limit 

Figure  3.4(a)  plots  the  proportional  limit  as  a  function  of  relative  density.  The 
magnitude  of  the  proportional  limit  ranges  from  1  MPa  to  3.1  MPa  over  all  the  testing 
conditions  and  has  no  consistent  trend  as  a  function  of  relative  density.  A  plot  of  the 
proportional  limit  as  a  function  of  confining  pressure  for  all  tests  conducted  at  moderate 
strain  rate  clearly  indicates  a  decrease  with  increasing  confining  pressure  [Fig.  3.4(b)]. 
This  trend  is  consistent  for  all  strain  rates  investigated.  Figure  3.4(c)  shows  the  effect  of 
strain  rate  on  the  proportional  limit  for  tests  conducted  at  ac~  0.1  and  10  MPa.  There  is 
a  tendency  for  the  proportional  limit  to  increase  with  strain  rate. 

3.3.4  Yield  Offset  at  10~^  Strain 

Figure  3.5(a)  plots  the  deviator  stress  at  the  10 yield  offset  as  a  function  of 
relative  density.  Over  all  testing  conditions,  the  yield  point  ranges  from  2.7  MPa  to  5.2 
MPa  with  no  apparent  trend  as  a  function  of  relative  density.  There  is  a  dear  trend  for  the 
offset  stress  to  decrease  with  increasing  confining  pressure  [Fig.  3.5(b)].  This  trend  is 
uniform  for  all  strain  rates  investigated.  As  a  fimction  of  strain  rate,  there  is  a  definite 
trend  for  the  yield  o&et  to  increase  with  strain  rate  with  a  power  law  coeffident  of 
approximatdy  11  [Fig.  3.5(c)]. 
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3.4  LARGE  STRAIN  BEHAVIOR 

This  section  discusses  the  large  strain  behavior  of  frozen  Manchester  Fine  Sand  as 
described  by  the  upper  yield  stress  and  corresponding  strain,  the  rate  of  post  upper  yield 
strain  hardening  or  softening,  the  peak  deviator  stress  and  corresponding  strain,  and  the 
volumetric  behavior. 

3.4.1  Overview  of  Large  Strain  Behavior 

The  influence  of  the  testing  variables  (Dr,  ffc  and  e)  on  the  nature  of  the  large 
strain  behavior  may  be  summarized  as  follows. 

1)  Change  in  Rdative  Density  (Dr)  at  Moderate  Strain  Rate 

As  illustrated  in  Fig.  3.6  (a)  and  (b),  increasing  Dr  causes  a  large  increase  in  the 
rate  of  post  upper  yield  strain  hardening  and  hence  peak  strength.  The  volumetric 
expansion  also  increases.  Note  that  the  deviator  stress  has  been  normalized  to  Quy 
since  Dr  does  not  affect  the  upper  yield  stress  (see  Section  3.4.2). 

2)  Change  in  Confining  Pressure  (^c)  at  Moderate  Strain  Rate 

Figures  3.7  (a)  and  (b)  show  that  increasing  ffc  also  causes  a  general  increase  in 
the  amount  of  strain  hardening  and  hence  peak  strength,  with  the  effect  being  more 
pronounced  at  high  densities.  The  amount  of  dilation  is  greatly  suppressed  with 
higher  levels  of  confinement.  Note  that  in  these  figures  the  deviator  stress  has  also 
been  normalized  since  ffc  has  a  minor  effect  on  Quy. 

3)  Change  in  Strain  Rate  (e)  at  varying  Density  and  Confinement 

Increasing  e  causes  a  large  increase  in  the  upper  yield  stress,  the  peak  strength 
and  increased  in  the  amount  of  dilation  [Fig.  3.8  (a),  (b)  and  (c)].  The  fastest  strain 
rate  causes  the  least  amount  of  post  yield  strain  hardening  or  even  strain  softening. 

3.4.2  Upper  Yield  Behavior 

The  upper  yield  region  represents  the  knee  of  the  stress-strain  curve  where  very 
significant  amounts  of  plastic  deformation  begin  to  occur.  This  region  is  emphasized  in 
these  data  because  the  use  of  accurate  on-spedmen  axial  strain  measurements  and  the 
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wide  lange  of  testing  variables  provide  the  first  comprehensive  evaluation  of  the  upper 
yield  behavior  for  firozen  soil.  Figure  3.9  illustrates  the  nature  of  the  stress-strain  data 
(plotted  up  to  axial  strain  levds  of  about  one  percent)  that  enabled  reliable  determinations 
of  the  upper  yield  stress  (Quy)  and  corresponding  strain  (ey)  as  a  function  of:  (a)  relative 
density;  (b)  confining  pressure  and  (c)  strain  rate.  (Note:  Section  5.4.1  of  Andersen  (1991) 
describes  the  technique  used  to  obtain  values  of  Quy  corrected  to  a  standard  initial 
specimen  height  of  7.5  cm;  Table  A.2  reports  both  the  measured  and  corrected  values). 

Figure  3.10  shows  that  the  magnitude  of  the  upper  yield  stress  is  independent  of 
rdative  density  and  Fig.  3.11  indicates  that  increasing  confining  pressure  causes  a 
consistent,  but  relatively  small,  decrease  in  Quy.  However,  increasing  strain  rate  causes  a 
very  significant  increase  in  upper  yield  stress  as  shown  by  the  results  in  Fig.  3.10  (Quy 
versus  rdative  density  for  all  tests)  and  in  Fig.  3.12  (Quy  versus  strain  rate  for  tests  with 
ffc  =  0.1  and  10  MPa).  The  power  law  coeffident  was  4.6  to  4.7  for  these  strain  rates.  The 
influence  of  testing  variables  on  the  magnitude  of  the  axial  strain  at  upper  yidd  follows  the 
same  basic  pattern:  that  is,  the  data  in  Fig.  3.13(a)  show  no  effect  of  rdative  density, 
whereas  the  data  in  Fig.  3.13(b)  show  significantly  larger  strains  with  increasing  strain 
rate.  In  fact,  the  magnitude  of  Cy  varies  almost  linearly  with  the  value  of  Quy 

Figure  3.14  compares  the  strain  rate  sensitivity  of  frozen  MFS  vrith  that  of 
isotropic  fine-grained  (0.7mm)  granular  ice  t^ted  at  —  7  1*C  (Hawkes  and  Mellor  1972) 

over  roughly  the  same  range  in  strain  rates.  Although  the  upper  yidd  stress  of  the  frozen 
sand  is  about  double  the  strength  of  ice,  both  systems  have  the  same  power  law  coeffident 
(n  =  4.7)  defined  by  the  slope  of  the  line  in  the  log  e  —  log  Q  plot.  This  significant 
observation  suggests  that  the  physical  mechanisms  controlling  the  yidd  behavior  of  frozen 
sand  (at  least  for  MFS)  might  be  similar  to  those  controlling  the  strength  of  polycrystalline 
ice.  In  any  case,  the  "strength"  of  the  sand  skeleton  must  not  contribute  significantly  to 
the  upper  yidd  stress  since  sand  density  did  not  influence  the  upper  yidd  behavior  and 
confining  pressure  had  a  minor  effect.  Further  research  is  warranted  to  determine  if 
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approaches  used  to  predict  the  behavior  of  composite  materials,  which  were  successful 
regarding  initial  modulus,  might  be  extended  into  the  upper  yield  region. 

3.4.3  Peak  Strength  Behavior 

Since  both  relative  density  and  confining  pressure  affect  the  post  upper  yidd 
behavior  of  frozen  MFS,  the  "strength"  of  the  sand  skdeton  becomes  important  at  large 
strains. 

Effect  of  Rel&tive  Density 

The  results  in  Fig.  3.15(a)  show  that  the  peak  strength  (Qp)  increases  linearly 
with  rdative  density  for  shearing  at  the  moderate  strain  rate.  The  data  also  show  that  the 
slope  increases  with  confining  pressure.  Since  the  upper  yidd  stress  remains  constant 
(Quy  =  3.1  ±  0.5  MPa  as  per  Fig.  3.10),  the  strength  increase  is  due  to  higher  rates  of  post 
yidd  strain  hardening  (e.g..  Fig.  3.6).  Increasing  density  produces  a  similar  trend  at  slow 
strain  rates  as  seen  in  Fig.  3.15(b).  Most  prior  data  also  show  a  linear  Qp  vs.  Dr 
rdationship,  e.g.,  Goughnour  and  Andersland  (1968);  Jones  and  Parameswaran  (1983). 
But  if  shearing  at  low  density  has  a  Type  I  curve,  such  as  test  FRS36  in  Fig.  3.8(a), 
increasing  density  may  only  lead  to  less  post  yidd  strain  softening.  Such  behavior  explains 
the  bilinear  Qp  vs.  Dr  rdationship  reported  by  Baker  and  Kurfurst  (1985). 

Figures  3.16(a)  and  (b)  plot  the  axial  strain  (cp)  at  the  peak  strength  versus 
rdative  density  for  tests  run  at  the  moderate  strain  rate  and  for  tests  run  at  the  slow  and 
fast  strain  rates,  respectivdy.  Increases  in  rdative  density  dther  have  little  effect  on  6p 
(generally  true  at  low  levds  of  confinement)  or  cause  a  rdatively  modest  decrease  in  Cp 
(always  true  at  high  confining  pressures). 

Effect  of  Confining  Pressure 

Figure  3.17  plots  peak  strength  versus  confining  pressure  at  three  rdative 
densities.  This  plot  was  devdoped  from  linear  regression  on  the  moderate  rate  test  data  in 
Fig.  3.15(a).  As  shown  in  the  figure,  strengthening  due  to  confinement  is  less  for  loose  sand 
than  for  dense  sand  and  tends  to  decrease  at  higher  stresses.  This  strengthening  is  always 
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accompanied  by  an  increase  in  the  axial  strain  at  peak  strength  (Fig.  3.16).  As  with 
rdative  density,  confinement  increases  the  amount  of  post  yield  strain  hardening,  especially 
for  dense  sand  (Fig.  3.7).  The  strengthening  efiect  is  similar  for  shearing  at  the  slow  rate 
[Fig.  3.15(b)].  However,  confinement  did  not  increase  the  strength  for  fast  shearing  of  loose 
sand  since  it  is  controlled  by  the  pressure  independent  upper  yield  stress;  e.g.,  tests  FRS36 
and  FRS40  in  Figs.  3.8(a)  and  3.8(c).  Hence  the  total  stress  friction  angle  of  frozen  MFS 
ranges  from  zero  (when  peak  =  upper  yidd  stress)  to  about  15*  (for  failure  at  large 
strains).  (Note:  sin  (|>  =  S/(2  +  S)  where  S  =  dQp/d^c)* 

Data  on  other  sands  also  show  large  variations  in  ()>  at  similar  confining  stresses. 
The  friction  angle  of  some  dense  frozen  sands  even  exceeds  20  to  25* ,  e.g.,  Chamberlain  et 
al.  (1972),  Alkire  and  Andersland  (1973)  and  Shibata  et  al.  (1985).  Thus  (j>  of  frozen  sand 
can  be  much  larger  than  measured  for  granular  ice  (say  ()>  =  10  ±  5*  at  moderate  to  fast 
strain  rates)  or  from  consolidated— undrained  tests  on  unfrozen  sand  (()>=::  5  to  10*). 

Effect  of  Strain  Rate 

Figure  3.18  summarizes  the  effect  of  strain  rate  on  the  peak  strength  of  loose 
and  dense  MFS  at  low  and  high  confining  pressures.  This  log-log  plot  presents  the  most 
complete  picture  of  the  influence  of  relative  density  and  confinement  on  the  rate  sensitivity 
of  a  frozen  soil.  The  corresponding  power  law  coefficients  (n)  that  denote  the  change  in  log 
c  per  change  in  log  Qp  range  from  about  6  to  17  (Note:  increases  in  n  reflect  decreases  in 
the  rate  sensitivity.)  Although  these  values  compare  well  vrith  prior  results  at  about 
— 10*C,  this  program  did  not  measure  a  rate  insensitive  region  (i.e.,  n  =  x)  such  as  observed 
by  Bragg  and  Andersland  (1980)  (for  unconfined  compression  tests  at  6  >  lO'Vsec  on 
medium-dense  sand)  and  by  Yuanlin  et  al.  (1988)  (for  unconfined  compression  tests  at  e  > 
3  X  lO'Vsec  on  dense  sand). 

Figure  3.18  shows  a  consistent  trend  for  the  power  law  coefficient  to  increase 
with  both  increasing  relative  density  and  increasing  confining  stress.  The  low  density-low 
confinement  condition  has  the  highest  rate  sensitivity  (n  =  6.1),  which  begins  to  approach 
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that  for  granular  ice  for  the  same  range  of  strain  rates  (e.g.,  n  =  4.7  from  Fig.  3.14).  From 
the  stress-strain  curves  in  Fig.  3.8(a),  one  sees  that  the  peak  strength  of  these  tests  either 
equals  or  only  modestly  exceeds  the  upper  yield  stress.  In  turn,  Section  3.4.2  suggests  that 
the  physical  mechanisms  controlling  the  yield  stress  of  frozen  MFS  might  be  similar  to 
those  controlling  the  strength  of  polycrystalline  ice. 

The  high  density— high  confinement  condition  has  the  lowest  rate  sensitivity 
(n  =:  16.7)  and  Fig.  3.8(b)  shows  that  the  peak  strength  for  these  tests  occurs  at  stress  and 
strain  levels  much  larger  than  those  at  the  upper  yidd  stress.  The  very  significant 
post— yield  strain  hardening  probably  reflects  a  substantial  contribution  from  the  frictional 
strength  of  the  sand  skeleton,  which  should  be  strain  rate  insensitive.  If  the  upper  yield 
stress  inherits  its  rate  sensitivity  from  the  pore  ice,  and  the  gain  in  strength  above  the 
upper  yield  point  is  much  less  rate  sensitive  because  of  sand  skeleton  friction,  it  follows 
logically  that  the  peak  strength  of  frozen  dense  sand  with  high  confinement  should  exhibit 
less  strain  rate  dependence. 

3.4.4  Volumetric  Behavior 

The  volumetric  strain  data  in  Figs.  3.6,  3.7  and  3.8  show  some  aspects  that  are 
consistent  with  the  drained  shearing  of  unfrozen  sand,  i.e.,  less  dilation  (expansion)  with 
lower  density  and  higher  confinement.  But  the  overall  behavior  is  very  different  from 
unfrozen  sand  in  several  respects: 

1)  High  confinement  does  not  cause  significant  contraction,  since  this  is  prevented  by 
the  ice  matrix  in  the  pores. 

2)  The  peak  strength  does  not  occur  at  the  maximum  rate  of  dilation,  dcy/de^. 

3)  The  amount  of  dilation  at  large  strains  is  much  larger  than  can  be  explained  by 
behavior  of  the  unfrozen  sand  skeleton.  In  other  words,  the  ice  matrix  causes 
additional  expansion.  Continuous  dilation  in  some  tests  at  low  density  and  low 
confinement  caused  final  densities  less  than  the  minimum  density,  i.e..  Dr  <  0%. 
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3.4.5  Summary  nf  Bffertu  nf  Rriativp  Demritv.  Confining  Pre8siire  and  Strain 
Rate 

The  test  program  on  MFS  has  genaated  the  most  complete  study  of  the  effects 
of  relative  density,  confining  pressure  and  strain  rate  on  the  strength-deformation 
properties  of  a  firozen  soil  that  is  available  in  the  literature.  Accurate  measurements  of 
behavior  from  very  small  to  very  large  strains  at  -9.5*  C  has  produced  some  unique  results, 
espedally  regarding  definition  of  the  upper  yidd  stress,  Quy.  The  magnitude  of  Quy  is 
independent  of  sand  density  (Fig.  3.10)  and  decreases  only  slightly  with  confining  pressure 
(Fig.  3.11).  This  suggests  that  the  upper  yidd  stress  is  not  affected  by  the  fiictional 
character  of  the  sand  skdeton.  Although  Quy  is  about  double  the  unconfined  peak  strength 
of  fine-grained  polycrystalline  ice,  both  systems  exhibit  similar  sensitivities  to  strain  rate 
(Fig.  3.14).  This  suggests  that  the  physical  mechanisms  controlling  the  upper  yidd 
behavior  of  frozen  sand  might  be  similar  to  those  controlling  the  strength  of  ice.  Section 
3.3.2  shows  that  the  initial  modulus  of  frozen  MFS  can  be  reasonably  predicted  by  treating 
the  system  as  a  composite  material,  i.e.,  "stiff"  sand  grains  embedded  in  a  "soft”  ice 
matrix. 

In  contrast  to  the  above  findings,  the  stress-strain  behavior  of  frozen  MFS  at 
large  strains  is  greatly  affected  by  the  fiictional  resistance  of  the  soil  skdeton,  except  for 
fast  shearing  of  loose  MFS  at  low  confinement,  where  failure  occurs  at  a  low  strain  and  the 
peak  strength  equals  the  upper  yidd  stress  [i.e.,  test  FRS36  in  Fig.  3.8(a)].  All  other 
conditions  (lower  e,  higher  Dr  and/or  higher  ctc)  lead  to  post  yidd  strain  hardening  and 
failure  at  larger  strains  (Fig.  3.2).  Under  these  conditions,  the  peak  strength  was  shown  to: 

1)  Increase  linearly  with  rdative  density,  with  the  rate  bdng  larger  with  confinement 
(Fig.  3.15). 

2)  Increase  nonlinearly  with  confinement,  with  the  rate  bdng  larger  at  higher  density 
(Fig.  3.17). 
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These  strengthening  effects  are  less  strain  rate  sensitive  than  granular  ice  (Fig.  3.18)  and 
presumably  reflect  a  greater  contribution  of  the  frictional  resistance  of  the  sand  skdeton  to 
the  peak  strength. 

3.5  EFFECT  OF  TEMPERATURE  ON  THE  BEHAVIOR  OF  FROZEN 

MANCHESTER  FINE  SAND 

All  but  one  of  Andersen’s  (1991)  tests  in  the  "prindpal"  program  were  conducted  at 
one  temperature  (— 9.5*C).  Since  Ashby  and  Duval  (1985)  showed  that  the  "normalized" 
primary  creep  behavior  of  ice  (within  the  ductile  deformation  regime)  was  the  same  for 
changes  in  applied  stress  (equivalent  to  changes  in  e)  and  changes  in  temperature,  the 
principal  program  focused  on  changes  in  strain  rate.  However,  at  the  suggestion  of  David 
Cole,  tests  at  different  temperatures  were  started  last  summer  (1991)  after  first  improving 
our  temperature  control  and  measurement  capabilities.  The  program  to  date  has  been 
restricted  to  six  moderate  strain  rate  tests  on  loose  sand  at  low  confinement  for  comparison 
with  the  effects  of  changing  strain  rate. 

Figure  3.19  plots  deviator  stress  versus  axial  strain  for  the  "new"  tests  (FRS  75  to 
80)  having  temperatures  ranging  from  -5.3*C  to  — 28*C,  plus  FRS44  as  being 
representative  of  tests  sheared  at  -9.5*  C.  Five  of  the  new  tests  exhibited  erratic  behavior 
at  large  strains  for  unknown  reasons;  the  dotted  curves  denote  these  questionable  data.  In 
addition,  FRS80  sheared  at  -28*  C  appears  to  have  failed  prematurdy,  again  for 
unexplained  reasons. 

The  expanded  scale  plot  in  Fig.  3.19(b)  indicates  that  temperature  has  very  little 
effect  on  small  strain  behavior.  Moreover,  the  values  of  Young’s  modulus  presented  in 
Table  A.ll  for  FRS75  to  80  give  E  =:  26.2  ±  3.8  MPa,  compared  to  26.6  ^  4.5  MPa  quoted 
in  Table  A.4  firom  the  principal  program.  This  dose  agreement  should  be  expected  since: 
Section  3.3.1  shows  that  modulus  is  independent  of  strain  rate;  and  Sinha  (1989)  shows 
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that  the  modaliu  of  graniilar  ice  increases  by  only  about  3%  as  the  temperature  decreases 
fiKMn-5  to-30*C. 

As  expected,  Fig.  3.19(a)  shows  a  significant  increase  in  both  the  upper  yidd  stress 
and  the  peak  strength  with  decreasing  temperature.  There  is  also  a  change  in  the  type  of 
curve  from  type  ni  for  shearing  at  T  =  —10  ^  5*C  to  t3rpe8  II  or  I  for  shearing  at  lower 
temperatures.  This  trend  is  very  similar  to  that  observed  for  the  influence  of  increasing 
strain  rate  on  the  behavior  of  frozen  MFS  at  low  density  and  low  confining  pressure  as 
shown  in  Fig.  3.8(a). 

Figure  3.20  plots  the  (corrected)  upper  yidd  stress  (Quy)  ^d  the  peak  strength 
(Qp)  as  a  function  of  temperature  for:  loose  sand  at  low  confinement  from  the  six  new 
tests  and  the  mean  values  from  three  tests  run  at  -9.5*  C;  and  dense  sand  at  high 
confinement  from  tests  in  the  prindpal  program.  It  should  be  emphasized  that  the 
temperature  in  test  FRS50  was  not  controlled  and  that  T  =  — 15*C  represents  an 
approximate  value.  The  low  Dr  —  low  ac  tests  show  an  approximatdy  linear  inaease  in 
both  Quy  and  Qp  with  decreasing  temperature  (exduding  test  at  — 28*C),  with  slopes  of  0.7 
and  0.6,  respectivdy.  A  summary  of  test  results  from  previous  programs  (Bourbonnais  and 
Ladanyi  1985)  also  shows  an  approximatdy  linear  rdationship  of  similar  slope  (dQ/dT  = 
0.7  ±  0.3)  for  the  peak  strength  of  various  frozen  sands  sheared  at  c  =  3  ^  2  x  lO'^/sec  over 
the  same  temperature  range.  The  "low"  value  of  Quy  for  the  test  FRS50  at  -15*  C  may 
reflect  uncertainty  in  the  actual  temperature  because  the  results  in  Section  3.4.2  suggest 
that  density  and  confinement  should  not  affect  dQuy/dT  (i.e.,  these  variables  did  not  affect 
the  rate  of  increase  in  Quy  with  increasing  strain  rate).  The  smaller  dQ/dT  slope  for  the 
peak  strength  at  high  Dr  —  high  Cc  needs  verification. 

As  noted  earlier  in  Section  3.5,  one  objective  of  testing  at  different  temperatures  is 
to  compare  it’s  effects  with  those  observed  for  changes  in  the  strain  rate.  An  important 
part  of  this  comparison  should  indude  an  evaluation  of  how  these  two  variables  influence 
the  behavior  of  frozen  MFS  vis-a-vis  that  of  granular  ice.  In  Section  3.4,  the  rate 
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sensitivity  of  frozen  MFS  was  directly  compared  to  the  rate  sensitivity  of  ice  using 
experimental  data  at  a  similar  temperature  (i.e.,  at  near  —10*  C).  However,  the  authors  did 
not  find  ice  data  to  enable  a  direct  comparison  of  its  temperature  sensitivity  to  that  of 
frozen  MFS  at  the  moderate  strain  rate  of  e  =  3  x  10'^/ sec.  For  example,  Jacka  (1984) 
presents  extensive  creep  data  at  temperatures  ranging  from  -fi*  C  to  —32.5*  C,  but  only  for 
strain  rates  slower  than  10~*  to  lO'Vsec,  i.e.,  well  within  the  ductile  (flow)  deformational 
r^me.  Consequently,  the  temperature  sensitivity  of  ice  was  estimated  using  the  relation 
between  strain  rate,  stress  and  temperature  presented  in  Glenn  (1955).  Dr.  Shyam  Sunder 
at  MIT  suggested  the  procedure  described  below  for  determining  the  "constants"  in  Glenn’s 
equations. 

Glen’s  (1955)  power  law  creep  equation  can  be  expressed  as  (e.g.,  Shyam  Sunder  and 
Wu,  1989) 

Qp  =  V  (e)^/“  . 3.2 

where  Qp  is  the  measured  peak  strength  (MPa),  i  is  the  strain  rate  (1/sec),  n  is  the  power 
law  coefficient  determined  from  the  experimental  data,  and  V  is  the  derived  temperature 
dependent  reference  stress  (MPa)  at  a  reference  strain  rate  of  unity.  In  turn,  the  Arrhenuis 
equation  was  used  to  predict  the  temperature  dependence  of  the  reference  stress 

V  =  Vo  exp(Q/nRT) . 3.3 

where  T  is  the  temperature  (*K),  Q  is  the  activation  energy  (assumed  equal  to  67  kJ/mol), 
R  is  the  universal  gas  constant  (8.31  J/mol  — *K),  and  Vq  is  a  temperature  independent 
constant  that  was  backcalculated  from  the  experimental  data. 


Three  sets  of  unconfined  compression,  constant  strain  rate  test  data  were  evaluated. 


Data  Set 

Reference 

Tcmp(*,g) 

e  Range 

Grain  Dia. 

A 

Hawkes  and 
MeUor  (1972) 

-7 

10-*  to  10-» 

0.7mm 

B 

Jones  (1982) 

-11.7 

1.4  X  10 -« to 

5.4  X  10-4 

1mm 

C 

Cole  (1987) 

-5 

10-5  to  10-3 

1.2mm 

Equation  3.2  was  first  used  to  determine  the  values  of  V  and  n  and  then  Eq.  3.3  was  used 
to  obtain  the  value  of  Vo,  with  the  following  results: 


Data  Set 

VfMPa^ 

& 

rl 

A 

-7 

39.7 

4.67 

60.5 

0.96 

B 

-11.7 

21.3 

6.48 

181 

0.95 

C 

-5 

60.9 

3.83 

23.4 

— 

The  above  values  of  n  and  Vo  were  used  with  Eq.  3.3  to  compute  the  values  of  V  at 
different  temperatures,  and  then  Eq.  3.2  to  calculate  Qp  at  e  =  3  x  lO'Vsec,  the  moderate 
strain  rate  used  for  the  frozen  MFS  tests.  It  is  recognized  that  this  approach  is 
approximate  since:  the  experimental  data  lie  within  the  ductile  to  brittle  transition  zone 
(hence  n  may  vary  with  temperature);  and  the  Arrhenius  equation  may  not  be  valid  at 
temper?  ^ures  above  — 10*C  (hence  Q  may  not  be  constant).  Hence  the  approach  was 
applied  to  three  data  sets  in  order  to  obtsun  a  range  in  the  predicted  temperature 
sensitivity  of  ice. 

Figure  3.21  compares  the  calculated  strength  of  ice  to  the  measured  upper  yield 
stress  (Quy)  peak  strength  (Qp)  of  frozen  MFS  at  various  temperatures.  The  figure 
plots  log  stress  versus  1/T  (*K)  since  this  format  produces  a  linear  relationship  for  ice 
based  on  Eqs.  3.2  and  3.3.  The  results  clearly  indicate  that  the  temperature  sensitivity  of 
frozen  MFS  is  very  different  from  that  predicted  for  ice:  the  log  stress  vs.  1/T  relationship 
for  frozen  MFS  is  highly  non-linear  and  has  a  much  larger  sensitivity  at  temperatures 
ranging  from  -6*  C  to  —15*  C. 

Thus  the  initial  "assumption”  that  the  effects  of  varying  temperature  and  of  varying 
strain  rate  on  the  behavior  of  frozen  sand  would  be  similar  is  certainly  not  true.  For 
example,  Fig.  3.14  showed  that  the  rate  sensitivity  of  Quy  for  frozen  MFS  and  of  Qp  for  ice 
(data  set  A)  were  the  same.  This  led  to  the  suggestion  in  Section  3.4.2  that  the  physical 
mechanisms  controlling  the  yield  behavior  of  frozen  MFS  might  be  similar  to  those 
controlling  the  strength  of  granular  ice.  That  hypothesis  needs  further  experimental 
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evaluation,  such  as  by  vaiying  the  strain  rate  at  different  temperatures.  Table  3.2  presents 
the  tentative  scope  of  further  studies  as  contained  in  the  Ladd  and  Germaine  (1991) 
proposal  to  ARO,  which  also  includes  testing  at  conditions  of  high  density  and  high 
confinement. 

3.6  EFFECT  OF  CONFINEMENT  DURING  FREEZING  ON  FROZEN  SAND 
BEHAVIOR 

The  1988  proposal  for  this  research  included  some  triaxial  tests  wherein  the 
confining  stress  would  be  applied  for  a  long  period  in  order  to  consolidate  the  frozen  soil 
and  thereby  increase  the  preshear  effective  stress  acting  on  the  sand  skeleton 
(freeze-consolidate  series  FC  in  Table  1.2).  There  are  two  problems  with  this  type  of 
testing: 

1)  It  would  require  very  long  times  (as  predicted  by  Ladanyi  1985)  and  thus  was  not 
feasible  with  only  one  high  pressure-low  temperature  triaxial  system; 

2)  It  does  not  simulate  in  situ  conditions  wherein  effective  stresses  act  on  the  soil 
skeleton  before  freezing  occurs. 

However,  the  concept  of  studying  the  influence  of  the  pre-freezing  effective  stress  on  frozen 
soil  behavior  remains  valid.  But  rather  than  first  freezing  and  then  consolidating  the  soil 
(FC  test),  the  process  will  be  reversed,  i.e.,  first  consolidate  and  then  freeze  the  soil 
(consolidate-freeze  =  CF  test).  As  noted  in  Section  4,  increasing  the  presbe&r  effective 
stress  should  cause  a  very  significant  increase  in  the  stiffness  of  the  soil  skeleton.  For 
example,  stress-strain  data  for  unfrozen  dense  MFS  in  Fig.  4.2  at  6a  =  2%  show  Q 
increasing  dramatically  from  0.7  to  5.5  MPa  as  a'c  increases  from  0.1  to  10  MPa.  This 
large  increase  in  sand  skeleton  stiffness  will  presumably  have  a  significant  affect  on  the 
strength-deformation  behavior  of  frozen  sand  at  strains  near  and  beyond  the  upper  yield 
r^on. 
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Our  review  of  the  literature  indicates  that  no  tests  have  been  run  to  compare  the 
behavior  of  frozen  8(^  as  measured  in  regular  "frozen— unconsolidated''  triaxial  compression 
tests  compared  to  tests  run  on  soil  that  was  frozen  after  application  of  confinement.  Singh 
et  al.  (1982)  did  subject  a  specimen  of  uniform  sand  (D50  =  0.4mm)  to  one-dimensional 
freezing  in  a  tiiaxial  cell  (consolidation  stress  =  0.56  kg/cm^  and  backpressure  = 
2.5kg/cm3)  with  measurements  of  the  change  in  pore  volume  during  freezing  and 
subsequent  thawing.  But  no  tests  were  run  on  frozen  sand  since  the  objective  was  to  show 
that  freezing— thawing  did  not  effect  the  cyclic  behavior  of  unfrozen  sand.  Goodman  (1975) 
measured  the  "viscoelastic  response"  and  yidd  stress  of  a  20—40  quartz  sand  that  was 
frozen  one-dimensionally  (T  =  — 7±rC)  after  first  applying  confining  pressures  of  125,  250 
and  500  psi  with  the  objective  of  simulating  deep  permafrost.  The  specimens  were  enclosed 
in  plasticized  PVC  tubing  during  application  of  the  confining  pressure  and  then  subjected 
to  axial  loading  creep  tests.  However,  since  the  yield  stress  was  shown  to  be  the  same  for 
both  frozen  and  unfrozen  sand,  the  writers  infer  that  the  loading  conditions  may  have  been 
closer  to  conSned  (1— D)  compression  than  triaxial  compression.  Sego  et  al.  (1982)  used  a 
triaxial  cell  for  one— dimensional  freezing  of  mortar  sand  under  a  consolidation  stress 
ff'c  =  0.08  MPa  in  order  to  study  the  influence  of  pore  fluid  salinity  on 
strength-deformation  behavior  at  varying  strain  rates  (T  =  — 7*C).  But  all  frozen 
specimens  were  transferred  to  another  cell  for  triaxial  compression  tests  at  essentially  zero 
confinement.  Thus,  the  program  did  not  provide  information  regarding  the  effect  of  the 
pre-freezing  effective  stress  on  the  behavior  of  frozen  sand. 

Test  results  from  the  current  program  may  indicate  how  the  pre— freezing  effective 
stress  affects  the  stress— strain  behavior  of  frozen  MFS.  Andersen  (1991)  compared  the 
stress— strain  behavior  of  dense  frozen  MFS  prepared  by  wet  tamping  (compaction)  rather 
than  by  multiple  sieve  pluviation  (MSP).  A  comparison  of  dense  (Dr  =  90%)  tests  under  2 
MPa  confinement  is  presented  in  Fig.  3.22.  This  figure  shows  that  compaction  caused  a 
significant  increase  in  the  rate  of  post  upper  yield  strain  hardening.  One  possible 
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explanation  for  this  behavior  is  that  the  compaction  process  "locked  in"  higher  effective 
stresses  and  that  these  stresses  were  preserved  during  freezing  and  hence  caused  higher 
preshear  effective  stresses  than  for  pluviated  specimens  (Note:  other  comparisons  at  ac  =  5 
and  10  MPa  also  showed  higher  post  upper  yield  sti&esses,  but  little  change  in  the  peak 
strength).  If  this  hypothesis  is  true,  even  rather  modest  changes  in  the  pre-ffeezing 
effective  stress  can  represent  an  extremdy  important  issue  about  which  no  information 
currently  exists. 


57 


Table  3.1  Scope  of  Testing  Prograa  on  the  Triaxial  Coq>re88ion  Behayior  of 
Frozen  Manchester  Fine  Sand 


RANGE 

IN 

RELATIVE  DENSITY 

(7.) 

CONFINING  PRESSURE  (MPa) 

0.1 

2 

5 

10  1 

Fast 

Mod. 

Slow 

Fast 

Mod. 

Slow 

Fast 

Mod. 

Slow 

Fast 

Mod. 

Slow 

20-40 

• 

• 

• 

• 

• 

• 

• 

. 

40-60 

• 

• 

• 

• 

■ 

60-  80 

• 

• 

• 

80- 100 

• 

• 

• 

• 

• 

• 

• 

•  1 

STRAIN  RATES  TESTED 

Fast  »  4  X  10‘</sec 
Moderate  »  3  x  10’ 5 /sec 
Slow  »  3  X  10’® /sec 
One  specimen  tested  at  1.2  x  10'® /sec 


TESTING  TEMPERATURE 

Principal  Program 
-9.55  ±  0.3’C 
New  Program  fto  date) 

Six  tests  at  temperatures  ranging 
from  -5  to  -28*C 

(Loose  sand  at  moderate  strain  rate  and  low  confinement) 


SPECIMEN  PREPARATION  TECHNIQUES 

Most  tests  prepared  by  Multiple  Sieve  Pluviation 
Seven  tests  prepared  by  Vet  Tamping 


SPECIMEN  END  CONDITIONS 

Most  tests  had  ice  cap  and  grease  ends  (low  friction) 
Six  tests  had  frozen  sand  and  grease  (higher  friction) 
One  test  had  emery  cloth 
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Tentative  Scope  of  Triaxial  Testing  on  IFS  at  Different  Teaperatnres 
("Conventional"  tests,  i.e.,  frozen  before  applying  confining  pressure) 


NORMALIZED  DEVIATOR  STRESS  NORMALIZED  OEVIATOR  STRESS 


a.  Low  Confining  Pressure 


AXIAL  STRAIN 


b.  High  Confining  Pressure 


AXIAL  STRAIN 

Note:  Curves  have  been  normalized  by  the  magnitude  of  the  upper  yield  stress 


Fignxe  3.2:  Overview  of  Tjrpes  of  Stress-Strain  Curves 
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Proportional  Limit  (MPa) 
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Figure  3.4(a):  Proportional  Limit  Versus  Relative  Density  for  Frozen 
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Figure  3.4(c):  Effect  of  Strain  Rate  on  Proportional  Limit  for  Frozen 

MFS 
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Figoie  3.5(b):  Effect  of  Confining  Pressure  on  Yield  Offset  Stress  for 

Frozen  MFS 


Figure  3.6(u): 


Normalized  Stress-Strain  Curves  Showing  the  Effect  of 
Relative  Density  for  Moderate  Strain  Rate  and  Low 
Confining  Pressure  for  Frozen  MFS 


VOLUMETRIC  strain  normalized  OEVIATOR  STRESS 


•  nr'.'  ■■  '.IK-.  ■£  ci  *1 


-  "  t.«i 


FRS19 
Dr  =  49.4% 


'e  =  Moderate 
ffc  =  10  MPa 


-0.5*  C 


liTi 


!  Effect  of 
and  High 


VOLUMETRIC  STRAIN  NORMALIZED  OEVIATOR  STRESS 


AXIAL  strain  '£  —  —9,5*  C 


Figtxre  3.7(a):  Normalized  Stress— Strain  Carves  Showing  the  Effect  of 

Confining  Pressure  for  Moderate  Strain  Rate  and  Low 
Native  Density  for  Frozen  MFS 
71 


VOLUMETRIC  I  RAIN  NORMALIZED  OEVIATOR  STRESS 


VOLUMETRIC  strain  OEVIATOR  STRESS  (MPa) 


AXIAL  STRAIN 


T  =  -9.5T 


Figuie  3.8(b): 


Stress-Strain  Carres  Showing  the  Effect  of  Strain  Rate 
for  High  Relative  Density  and  High  Confining  Pressure 
on  Frozen  MFS 


VOLUMETRIC  STRAIN  DEVIATOR  STRESS  (MPo) 
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Figure  3.8(c): 


Stress— Strain  Corves  Showing  the  Effect  of  Strain  Rate 
for  Low  Relative  Density  and  High  Confining  Pressure 
on  Frozen  MFS 


Figure  3.9(a):  Stress— Strain  Curves  Showing  Upper  Yield  Behavior  of 

Frozen  MFS  at  Varying  Relative  Density 


FRS53 
:  5.0  MPa 


Hgnre  3.9(b):  Stress-Strain  Curves  Showing  Upper  Yield  Behavior  of 

Frozen  MFS  at  Varying  Confining  Pressure 
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Flgme  3.9(c):  Stress-Strain  Curves  Showing  Upper  Yield  Behavior  of 

Frozen  MFS  at  Varying  Strain  Rate 
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Figure  3.10:  Corrected  Upper  Yield  Stress  Versus  Relative  Density 

for  Frozen  MFS 
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Fignie  3.15(a):  Peak  Deviator  Stress  Versus  Relative  Density  at 

Moderate  Strain  Rate  for  Frozen  MFS 
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Figure  3.15(b):  Peak  Deviator  Stress  Versus  Relative  Density  for  Range 

in  Strain  Rates  and  Confining  Pressures  for  Frozen 


Figme  3.16(a):  Axial  Strain  at  Peak  Deviator  Stress  Versus  Relative 

Density  for  Moderate  Strain  Rate  and  Range  in 
Confining  Pressures  for  Frozen  MFS 
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Figue  3.16(b):  Axial  Strain  at  Peak  Deviator  Stress  Versus  Relative 

Density  for  Range  in  Strain  Rates  and  Confining 
Pressures  for  Frozen  MFS 


Peak  Deviator  Stress  (MPa) 


Li-  CO 


90 


Strength  (Q 


Flgme  3.21:  Comparison  Between  Measured  Upper  Yield  Stress  and 

Peak  Strength  of  Frozen  MFS  and  Calculated  Peak 
Strength  of  Ice  at  Moderate  Strain  Rate 


VOLUMETRIC  STRaiU  DEVIATOR  STRESS.  MPo 


T  =  -9.5’C 

AXIAt  STRAiri 


Stress— Strain  Curves  to  Show  the  Effect  of  Specimen 
Preparation  Technique  on  Frozen  MFS 


Figure  3.22; 
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4.  BEHAVIOR  OF  UNFROZEN  MANCHESTER  FINE 
SAND  IN  TRIAXIAL  COMPRESSION 


4.1  SCOPE  OF  TESTING 

The  unfrozen  Manchester  Fine  Sand  (MFS)  testing  program  consisted  of 
isotropically  (hydrostatically)  consolidated— undrained  (CIUG)  and  consolidated-drained 
(CIDC)  triaxial  compression  tests.  The  majority  of  the  tests  (80%)  were  undrained.  The 
program  was  performed  in  two  series,  Series  A  and  Series  B,  as  summarized  in  Table  4.1 
and  discussed  below.  Of  the  total  56  tests  set-up,  only  38  (68%)  provide  useful 
information  for  analysis.  The  reasons  for  the  deletion  of  certain  tests  are  also  described. 
Spedmen  history  and  stress-strain  data  from  the  successful  tests  are  tabulated  in 
Appendix  B. 

4.1.1  Series  A  Tests 

For  Series  A,  16  tests  were  set-up  with  a  range  in  effective  confining  stresses 
(a'c)  of  0.1,  1,  2,  4  and  5  MPa.  All  except  one  test  were  sheared  under  undrained 
conditions.  Preshear  relative  densities  (Dr)  ranged  from  14  to  106  percent.  A  nominal 
strain  rate  of  2.3  x  lO'Vsecond  was  used  for  all  tests. 

Of  the  16  tests  set-up,  only  six  were  successful.  Several  factors  contributed  to  this 
low  success  rate.  Some  of  the  difficulties  involved  developing  new  experimental  techniques, 
e.g.,  lubricated  end  platens,  on-spedmen  axial  strain  measurements,  internal  axial  load 
cell,  high  capadty  pressure/volume  control  units,  and  spedmen  preparation  techniques. 
The  spedmen  preparation  techniques  for  Series  A  tests  induded  open-air  raining  of  MFS 
into  the  mold  (i.e.,  the  spedmen  former  as  illustrated  in  Fig.  2.10),  vibratory  densification 
of  rained  spedmens  and  "undercompaction”  (a  "wet  tamping"  method  where  the  spedmen 
is  formed  by  compacting  layers  inside  the  mold,  see  R.S.  Ladd  1974).  None  of  these  tests 
were  prepared  using  multiple  sieve  pluviation.  These  preparation  techniques  created 
non-unformities  in  the  test  spedmens  that  lead  to  numerous  erratic  and  questionable  test 
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results.  The  most  common  anomalies  in  the  testing  results  fell  into  two  areas:  1)  the  axial 
load  on  the  specimen  would  suddenly  drop  to  near  zero,  which  would  coincide  with  a 
sudden  jump  in  measured  pore  pressure;  and  2)  a  predominant  failure  surface  would 
devdop  prematurdy  in  the  specimen  during  shear,  causing  the  measured  pore  pressures 
and  axial  stress  to  vary  unpredictably  during  shear. 

Additional  problems  encountered  during  Series  A  tests  included: 

1)  Membrane  breakage  and/or  leakage  during  consolidation  or  shearing; 

2)  Malfunction  of  the  on-spedmen  ACDTs  thereby  losing  small  strain  information; 

3)  Indeterminate  volume  changes  during  the  (attempted)  simultaneous  increase  of  the 
cell  and  back  pressure  for  overnight  saturation.  In  some  cases,  the  set— up  effective 
stress  was  lost  during  this  manually-controlled,  pressure— up  procedure. 

Because  of  these  problems,  many  spedmens  rarely  reached  the  large  strain  region,  i.e.,  15 
to  25%  axial  strain. 

4.1.2  Series  B  Tests 

Forty  tests  were  performed  in  the  second  set  of  tests.  Series  B.  All  spedmens 
for  Series  B  were  prepared  using  multiple  sieve  pluviation.  It  is  believed  that  this  one 
change  in  testing  procedure  was  responsible  for  dramatically  increasing  the  success  rate  for 
this  series.  Of  the  40  tests  performed  (32  CIUC  and  8  CIDC),  32  tests  were  successful. 
Preshear  relative  densities  varied  from  44  to  104%  and  effective  confining  pressures  were 
0.1,  2,  5  and  10  MPa,  except  for  a  test  at  7.5  MPa  and  another  at  11  MPa.  Strain  rates 
varied  from  2.3  to  3.5  x  lO'^/sec.  The  Series  B  tests  not  only  had  a  higher  success  rate 
than  for  Series  A,  but  also  reached  higher  strains  (25%)  with  better  stability  than  those  in 
Series  A.  Hence,  the  analysis  that  follows  is  based  largely  on  the  results  from  the  Series  B 
tests. 
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4.2  UNDRAINED  STRESS-STRAIN  AND  EFFECTIVE  STRESS  PATH  BEHAVIOR 

The  majority  of  useful  tests  are  CIUC  tests  (32  of  the  38);  therefore,  the  discussion 
will  focus  on  the  undrained  behavior  of  MFS.  Figure  4.1  presents  a  plot  of  qm  =  0.5(ori  - 
ffi)  versus  p'm  =  0.5(a'  i  +  or'z),  where  qn  and  p'a  correspond  to  conditions  at  the  point  of 
maximum  obliquity,  Raaz  =  (Note:  Fig.  4.1  also  includes  data  from  six 

CIDC  tests.)  Most  of  the  undrained  tests  results  have  effective  friction  angles  within  a 
range  of  32  to  37  degrees.  The  failure  envelope  has  a  noticeable  curvature  (concave 
downward),  which  is  typical,  e.g.,  as  observed  by  Lee  and  Seed  (1967)  ffom  drained  triaxial 
compression  tests  on  the  Saaemento  River  Sand. 

Figure  4.2  plots  deviator  stress  and  excess  pore  pressure  versus  strain  for  four  CIUC 
tests  at  different  effective  confining  pressures  (0.1,  2,  5  and  10  MPa)  on  dense  specimens 
having  almost  identical  preshear  void  ratios  (e  =  0.606  ±  0.007,  Dr  =  91  ±  2%).  Test  22 
(a'c—  0.1  MPa)  exhibits  a  strong  dilative  response  as  indicated  by  the  development  of 
large  negative  pore  pressures  which  strengthen  the  sand  as  straining  continues.  In  contrast, 
test  23  {<T'c  =  10  MPa)  exhibits  a  strong  contractive  response,  i.e.,  development  of  large 
positive  pore  pressures.  This  behavior  is  typical  of  similiar  comparisons  at  higher  preshear 
void  ratios  (i.e.,  lower  densities).  Figure  4.3  presents  the  effective  stress  paths  of  these 
same  four  tests.  The  dilatancy— hardening  response  of  test  22  causes  its  effective  stress 
path  to  "climb  up"  the  failure  envelope  and  requires  almost  20%  strain  to  reach  its  peak 
strength.  As  the  effective  confining  pressure  increases,  the  stress  paths  move  more  to  the 
left  before  climbing  the  failure  envelope,  until  at  10  MPa  the  sand  exhibits  only  contractive 
behavior.  In  fact,  test  23  may  have  failed  prematurely  (this  specimen  had  significant 
sliding  at  the  top  cap  and  the  after-shear  form  was  S-shaped). 

Figure  4.4  shows  the  influence  of  preshear  density  on  stress-strain  behavior  for 
three  CIUC  tests  having  an  effective  confining  pressure  of  2  MPa.  As  expected,  the  looser 
the  specimen,  the  stronger  is  its  contractive  response  during  shear  (generation  of  higher 
excess  pore  pressures)  and  the  lower  the  undrained  strength.  This  trend  in  behavior  is 


more  readily  seen  &om  a  plot  of  the  effective  stress  paths  in  Fig.  4.5.  The  loosest  specimen 
(test  39)  has  the  most  initial  contraction  and  little  strain  hardening,  whereas  the  densest 
specimen  (test  7)  shows  the  opposite. 

4.3  BEHAVIOR  IN  TERMS  OP  STEADY  STATE  LINE  AND  STATE 

PARAMETER  « 

Been  and  Jefferies  (1985)  introduced  a  concept  that  helps  to  unify  the  effects  of 
density  (Dr)  and  confining  stress  {a'c)  the  behavior  of  sands.  The  state  parameter  $  is 
a  measure  of  the  physical  condition  (state)  of  a  sand  in  terms  of  its  initial  void  ratio  and 
initial  state  of  stress  with  respect  to  an  ultimate  or  steady  state  condition.  They 
proposed  the  use  of  the  mean  normal  stress  (i.e.,  the  first  invariant  (Ii  =  a'oct)  of  the 
stress  tensor)  as  a  suitable  measure  of  the  stress.  Ii  is  defined  as  i  +  a' 2+  3)- 

Been  and  Jefferies  (1985)  assume  that  all  specimens  of  a  given  sand  tend  to 
approach  the  same  steady  state  irrespective  of  the  initial  state.  They  adopt  the  concept  of 
"the  Steady  State  of  Deformation"  as  proposed  by  Poulos  (1981)  as  the  ultimate  condition 
of  the  sand.  They  explain,  however,  that  their  concept  of  the  state  parameter  9  "does  not 
depend  on  the  nature  of  the  sand  structure  at  the  steady  state.  Rather  it  depends  on  there 
being  a  unique,  repeatable  particle  arrangement  at  the  steady  state  condition". 

Figure  4.6  illustrates  how  the  state  parameter  is  defined.  The  ordinate  is  the  void 
ratio  and  the  abscissa  is  the  mean  normal  stress  (Ii)  on  a  log  scale.  The  '9  parameter  is 
defined  as  tne  vertical  distance  from  the  initial  state  to  the  Steady  State  Line  (SSL) 
expressed  in  units  of  void  ratio.  A  sand  which  has  an  initial  state  that  plots  above  the  SSL 
has  a  +9  and  one  that  plots  below  the  SSL  has  a  —’I'.  Been  and  Jefferies  (1985  and  1986) 
report  that  many  sands  with  the  same  ^  exhibit  the  same  general  behavior  (e.g.,  the  same 
peak  drained  0'  for  —9  values  and  the  same  peak  undrained  strength  ratio  for  +'1'  values) 
irregardless  of  the  absolute  magnitude  of  the  Dr  or  a'c.  Thus  9  can  apparently  unify  the 
effects  of  both  Dr  and  a'c- 
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To  define  the  Steady  State  Line,  Been  et  al.  (1991)  state  that  tests  must  achieve  a 
condition  of  constant  volume  (if  drained)  or  constant  pore  pressure  (if  undrained)  during 
continued  shear  deformations  at  constant  shear  stress.  Based  on  this  criterion  and 
concentrating  on  the  undrained  test  results,  it  was  judged  that  16  tests  either  reached  or 
were  near  the  SSL  during  shear.  Where  appropriate,  the  selected  steady  state  condition  is 
shown  on  the  stress-strain  curve  and  the  effective  stress  path.  Figure  4.7  plots  the  mean 
effective  stress  (Ii  =  ff'oct)  versus  void  ratio  for  these  tests.  Appendix  B  presents 
tabulated  results  of  the  conditions  at  steady  state  for  these  tests.  A  Steady  State  Line  with 
A  =  0.253  was  determined  using  a  linear  regression  analysis  on  this  set  of  data  (r^  =  0.95) 
and  is  also  shown  in  Fig.  4.7.  It  is  important  to  note  that  the  Steady  State  Line  was 
developed  from  tests  having  both  negative  (moving  from  left  of  the  line)  and  positive 
(moving  from  right  of  the  line)  values  as  shown  by  the  arrows  in  the  figure. 

The  state  parameter  can  be  correlated  to  undrained  shear  parameters  such  as  the  A 
parameter  at  the  peak  strength  (Af)  and  the  undrained  strength  ratio  (qt/a'c)  as  shown  in 
Fig.  4.8  and  4.9,  respectively.  These  indicate  an  excellent  correlation  for  both  Af  and 
qf/o"c. 

The  concept  and  use  of  the  SSL  and  'I'  parameter  has  experienced  strong  debate  in 
the  literature,  e.g.  see  Been  and  Jefferies  (1986).  The  writers  agree  that  the  choice  of 
exactly  where  the  steady  state  occurs  can  be  subject  to  error.  Whereas  test  23  in  Fig.  4.2 
appears  to  have  reached  a  steady  state  condition  (before  undergoing  strain  softening, 
probably  due  to  cap  slippage),  tests  7  and  15  in  that  figure  did  not  exhibit  continued 
straining  at  constant  shear  stress  (although  the  final  points  were  selected).  Likewise,  the 
final  points  were  selected  for  tests  11  and  39  in  Fig.  4.4  even  though  the  shear  stress  was 
still  changing. 

Other  factors  also  influence  evaluation  of  the  SSL.  These  include  determination 
(measurement)  of  the  correct  void  ratio  of  a  test,  specimen  non— uniformities  (especially 
true  for  frictional  end  conditions)  and  grain  crushing  during  shear  which  is  applicable  for 
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tests  at  higher  confining  or  shear  stresses.  In  addition,  a  large  amount  of  strain  is  usually 
required  to  obtain  a  steady  state  condition  in  MFS;  therefore,  specimen  stability  on 
lubricated  ends  becomes  a  critical  issue.  However,  in  spite  of  these  potential  problems,  the 
data  in  Fig.  4.7  are  reasonably  consistent,  and  the  selected  SSL  gives  an  excellent 
corrdation  for  Af  and  Qf/a'c  (Fig-  4.8  and  4.9).  Nevertheless,  some  additional  study  of  the 
undrained  behavior  of  MFS  is  warranted. 

4.4  DRAINED  STRESS-STRAIN  AND  VOLUBfETRIC  BEHAVIOR 

A  study  of  the  drained  behavior  of  MFS  was  more  recently  initiated.  Six  successful 
tests  were  performed,  all  at  constant  p'  =  0.5(a'i  -f  a' 3).  Figure  4.10  shows  preliminary 
stress  and  volumetric  strain  results  of  four  tests  at  different  effective  confining  stresses  but 
with  very  similiar  preshear  densities.  The  test  at  a'c  =  0.1  MPa  exhibits  the  most 
expansion  (dilation)  during  shear.  (Note:  problems  occurred  after  about  8%  axial  strain.) 
As  the  confinement  increases,  the  strength  increases  and  the  tendency  of  the  specimen  to 
dilate  decreases.  Figure  4.11,  which  shows  the  stress  paths  for  these  tests,  illustrates  the 
curvature  of  the  failure  envelope  with  increasing  stress. 

One  CIDC  test  (No.  37)  meets  Been  et  al.  (1991)  definition  of  steady  state.  The 
steady  state  condition  of  this  test  is  plotted  on  Fig.  4.12  along  with  the  SSL  as  defined  by 
the  CIUC  tests  on  MFS.  Although  not  conclusive,  this  CIDC  "state  point"  appears  to 
validate  the  use  of  one  SSL  for  both  drained  and  undrained  MFS  tests.  However,  more 
testing  is  necessary  before  a  more  definitive  conclusion  can  be  made. 


100 


Table  4.1  Scope  of  Testing  Prograa  (m  Thfrozen  landiester  Fine  Sand 

in  Tfiaxial  Co^ression 


SERIES  A 

15  Undrained  Tests  (GIUC)  -  6  successful 
1  Drained  Test  (CIDC)  -  0  successful 

SERIES  B 

32  Undrained  Tests  (CIUC)  -  26  successful 
8  Drained  Tests  (CIDC)  -  6  successful 


VARIABLES  OF  THE  SUCCESSFUL  TESTS 


RANGE  IN  STRAIN  RATES 
«  2.3  X  10‘*/sec  to  3.7  x  10'5/sec 

SPECIMEN  PREPARATION  TECHNIQUES 
Open  Air  Raining  1 

Vibratory  Densification  V  Series  A 

Under  Compaction  (Vet  tamping) 

Muliple  Sieve  Pluviation  -  Series  B 

SPECIMEN  END  CONDITIONS 
Silicone  lubricant  only 
Rubber  Membrane  only 
Rubber  membrane  and  silicone  lubricant 
Rubber  membrane  and  high  vaccuum  grease 
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figure 4.3:  Effect  of  ConfioingPressure on  Effective  Stress  Paths  of 

Unfrozen  Dense  MFS  for  CIUC  Tests 
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Fignie  4.5:  Effect  of  Density  on  Effective  Stress  Paths  for  CIUG 

Tests  on  Unfrozen  MFS  with  <t\  —  2.0  MPa 
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Figure  4.6: 


Definition  of  State  Parameter,  9  (after  Been  and 
JdEferies  1985) 


-  parameter  at  Failure 


Flgnie  4.8:  Af  Versus  4^  Parameter  from  CIUC  Tests  on  Unfrozen 


Undrained  Shear  Strength  Ratio 
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Figure  4.9:  Normalized  Undrained  Shear  Ratio  Versus  'i'  Parameter 

from  CIUC  Tests  on  Unfrozen  MFS 
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Figure  4.10: 


Stress— Strain  and  Volumetric  Strain  Response  for  CIDC 
Tests  on  Dense  Unfrozen  MFS 


TEST 


Figure  4.11:  Effective  Stress  Paths  of  CIDC  Test  Results  on  Dense 

Unfrozen  MFS 


Void  Ratio 


Figure  4.12:  Steady  State  Line  and  CIDC  Test  Result  for  Unfrozen 


5.  APPUCATION  OF  DILATANCY-HARDENING  MODEL 

5.1  LADANYI’S  MODEL 

Ladanyi  (1985)  developed  a  dilatancy— hardening  theory  to  predict  the  strength  of 
frozen  dense  sand  and  Ladanyi  and  Morel  (1990)  present  experimental  results  to  check  the 
validity  of  the  model.  The  basic  concept  assumes  that  the  effective  stress  path  followed  by 
the  soil  skeleton  in  a  frozen  sand  will  be  the  same  as  that  for  unfrozen  sand  provided  that: 

1)  both  systems  are  subjected  to  the  same  strain  path  (axial  loading  at  constant  confining 
pressure  and  constant  volume);  and  2)  the  sand  skeleton  in  both  systems  starts  &om  the 
same  "state"  (density  and  effective  stress).  Unffozen  dense  sand  sheared  undrained  at  low 
confinement  (i.e.,  — ’^)  tends  to  dilate  (expand  in  volume).  This  causes  development  of 
negative  pore  pressures  that  increase  the  effective  stresses  acting  on  the  sand  which  results 
in  increased  frictional  resistance.  The  dilatancy-hardening  model  makes  the  following 
assumptions  (Ladanyi  and  Morel  1990): 

1)  All  of  the  pore  water  in  the  sand  is  considered  to  be  frozen. 

2)  The  behavior  of  the  sand  is  a  function  of  the  state  parameter  $  as  introduced 
by  Been  and  Jefferies  (1985)  and  described  in  Section  4. 

3)  No  consolidation  occurs  prior  to  shear,  i.e.,  the  tests  are  considered  as 
"unconsolidated".  The  shear  starts  ffom  a  known  density  and  confining 
pressure. 

4)  So  long  as  the  pore  ice  is  continuous  and  unbroken  during  shear,  the  sand  will 
behave  as  "undrained".  When  the  pore  ice  fails,  i.e.,  breaks  up,  the  sand 
behaves  as  "drained".  This  breakup  of  the  ice  has  been  associated  with  axial 
strains  of  1  to  2%  in  the  literature. 

5)  As  long  as  the  test  behaves  undrained  (i.e.,  before  the  pore  ice  becomes  broken), 
the  dilatancy-hardening  principles  established  for  unfrozen  sand  are  also 
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applicable  to  frozen  sand  provided  that  proper  account  is  made  for  the  tensile 
strength  of  the  pore  ice  matrix. 

Before  further  description  of  the  dilatancy— hardening  model,  the  definition  of  the 
critical  confining  pressure,  a'zcr,  is  needed.  Seed  and  Lee  (1967)  and  others  have  observed 
that  if  a  dense  sand  is  sheared  in  an  undrained  mode  with  no  back  pressure  in  a  state  where 
it  exhibits  a  tendency  to  dilate,  i.e.,  it  has  a  large  — $  parameter,  the  pore  water  pressure 
decreases  and  the  effective  stresses  increase  on  the  sand  skeleton  until  one  of  two  conditions 
occurs.  Either  the  pore  water  pressure  reaches  the  "tensile  strength"  of  water,  i.e.,  it 
decreases  sufficiently  to  cause  cavitation,  or  the  effective  stresses  increase  until  the  "state 
point"  reaches  the  Critical  Void  Ratio  Line  (CVRL)  or  alternatively,  the  Steady  State  Line 
(SSL).  The  transition  confining  pressure  between  these  two  condition  is  termed  the  critical 
confining  pressure,  ff'  scr-  Seed  and  Lee  (1967)  define  the  critical  confining  pressure  as  that 
confining  pressure  for  a  given  preshear  void  ratio  which  results  in  no  net  volume  change  at 
the  peak  strength  in  a  drained  triaxial  compression  test  on  unfrozen  sand  (or  no  net  change 
in  pore  piessure  at  failure  for  an  undrained  test).  As  used  by  Ladanyi  et  al.,  a'zcr  is  the 
initial  confining  pressure  that  will  result  in  no  change  in  pore  pressure  at  failure  for 
undrained  shear  of  unfrozen  sand  at  a  given  initial  density  (or  void  ratio). 

Although  Ladanyi  et  al.  presented  stresses  in  the  form  Q  =  (ai  —  03)  and  croct  = 
i(^i  +  +  0-3),  its  application  will  be  illustrated  using  the  "MIT"  format:  q  =  0.5(oi  — 

ffz)  and  p  =  0.5(ai  +  03)  as  represented  by  Fig.  5.1.  Summarizing  this  figure,  undrained 
shear  of  unfrozen  sand  starting  from  P3  =  a'zcr  will  end  up  at  point  B.  For  frozen  sand 
starting  from  the  same  P3  =  o’'3cr)  Ladanyi  et  al.  add  qi  (the  shear  strength  of  the  ice 
matrix)  and  thereby  obtain  point  C.  To  explain  further,  consider  the  strength  of  the 
unfrozen  soil.  Refering  to  Fig.  5.1,  the  effective  stress  envelope  starts  at  the  origin  and  has 
a  slope  equal  to  Kf  =  (a'3/(r'i)at  failure-  '^he  cavitation  domain  is  represented  by  a  total 
stress  failure  envelope  with  an  initial  slope  identical  to  the  effective  stress  envdope,  but 
translated  to  the  left  by  an  amount  equal  to  the  tensile  strength  of  the  pore  fluid.  The 
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transition  point  between  the  cavitating  and  noncavitating  regimes  occurs  for  a  confining 
pressure  of 

<rc“®’^3cr  —  T . 5.1 

where  T  is  the  tensile  strength  of  the  pore  fluid  (T  =  Tw  for  water  =  1  atm  and  T  =  Ti  for 
ice  »  1  atm).  For  a  confining  pressure  greater  than  scr  —  T,  the  total  stress  envdope  is 
horizontal  with  a  value  of  q  =  B.  The  predicted  strength  of  the  soil  skeleton  in  this  region 
is  therefore  solely  controlled  by  the  magnitude  of  the  a'  scr  uid  the  effective  friction  angle 

Since  pore  ice  can  support  both  tensile  and  shear  stresses,  Ladanyi  proposed  that  the 
shear  strength  of  the  ice  be  added  directly  to  the  shear  strength  of  the  sand  skeleton. 
Thus,  the  total  stress  line  (TSL)  for  the  frozen  sand  plus  ice  (FS  +  ice)  is  o&et  vertically 
by  a  distance  equal  to  qi  above  the  TSL  (FS)  line. 

Pertinent  equations  for  predicting  the  shear  strength  (q^g)  of  "relatively  dense 
frozen  sand"  from  Ladanyi’s  model  are  as  follows. 

Noncavitation  Case 

In  this  case,  az  is  greater  than  or  equal  to  <r'  scr  —  Ti,  where  Ti  =  tensile  strength  of 
ice.  The  shear  strength  of  the  sand  skeleton,  qss,  will  equal  0.5  a'zcr  (Rf  —1),  where 
Rf  =  (a'i/(T'3)at  failure  (i  e.,  1/Kf).  Therefore,  the  shear  strength  of  frozen  soil,  q__,  can 
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be  derived  by  adding  the  shear  strength  of  the  sand  skeleton  to  the  shear  strength  of  the  ice 
matrix 

qj.g  —  0*5  ff'  3cr  (Rf  ~  1)  +  Qi . 5.2 

For  this  case,  where  the  initial  confining  pressure  is  greater  than  or  equal  to  Pi  in  Fig.  5.1, 
the  model  predicts  a  constant  strength;  i.e.,  q^g  does  not  increase  with  increasing 
confinement. 

Cavitation  Case 

Here,  az  is  less  than  a'  za  —  Ti;  therefore,  the  shear  strength  of  the  sand  skeleton  is 
defined  as  0.5  {az  +  Ti)(Rf  —  1).  The  shear  strength  of  frozen  soil  is  now  given  by 
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4p3  —  0.5(<rj  +  Ti)(Bf  —  1)  +  qi . 5.3 

U  the  initial  confining  piessnie  is  less  than  Fi,  the  modd  predicts  that  q_-  increases  with 

FS 

increasing  confinement.  The  slope  of  this  total  stress  envdope  equals  Kf,  the  effective 
stress  failure  envdope  of  the  sand  skdeton. 


5.2  APPUCATION  OF  LADANYI’S  MODEL  TO  FROZEN  MFS 

As  stated  by  Ladanyi  and  Mord  (1990),  "undrained  tests  with  pore  pressure 
measurements  would  have  been  more  appropriate"  than  thdr  use  of  consolidated-drained 
triaxial  compression  tests.  Section  4  presented  the  results  of  triaxial  compression  tests 
characterizing  the  CIUC  behavior  of  MFS  as  a  function  of  density  and  consolidation 
(confining)  pressure.  Also  presented  was  a  Steady  State  Line  (SSL)  for  MFS,  thus  allowing 
the  measured  undrained  shear  behavior  to  be  evaluated  in  terms  of  the  state  parameter 

Figure  4.8  plots  Af  [Au/(<ri  —  (7a)f]  versus  r/f  and  shows  that  if  equals  —  0.07  at 
Af  ss  0;  i.e,  where  there  is  no  change  in  pore  pressure  at  failure.  Given  this  result,  one  can 
define  a  line  representing  the  critical  confining  pressure  a*  acr  versus  preshear  void  ratio  by 
offsetting  the  SSL  (presented  in  Figure  4.7)  downward  0.07  along  the  void  ratio  axis.  This 
line  is  presented  in  Fig.  5.2.  As  illustrated  in  the  figure,  for  dense  frozen  MFS  (Dr  =  95%), 
this  gives  a'  acr  =  2.37  MPa. 

Figure  5.3  compares  predicted  strengths  from  Ladanyi’s  model  to  measured  upper 
yidd  stress  (Quy)  auid  peak  strength  (Qp)  data  on  frozen  MFS  at  confining  pressures  of  0.1 
and  10  MPa.  Assumed  ice  properties  for  the  predictions  were  Ti  =  2.0  to  2.1  MPa  (taken 
from  Hawkes  and  Meilor  1972)  and  Qi  equal  to  those  illustrated  on  Fig.  5.3  (taken  from 
Jones  1982).  Since  Ti  is  very  close  to  the  critical  pressure  the  strength  of  the  sand  skeleton 
was  computed  using  Eq.  5.3  (cavitation  case)  for  =  0.1  MPa  and  Eq.  5.2  (non  cavitation 
case)  for  =  10  MPa.  For  low  confinement  (oc  =  0.1  MPa)  and  tr'scr  =  2.37  MPa, 
Qm  =  4.74  and  5.65  for  =  32'  and  35* ,  respectively.  For  high  confinement 
(ffc  =  10  MPa),  Qss  =  5.34  and  6.38  MPa  for  =  32*  and  35* ,  respectively.  In  Fig.  5.3, 
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the  hatched  zone  is  the  calculated  Q-.  for  (Tc  =  0.1  MPa  with  ice  strengths  taken  from 
nnconfined  tests  (ac  =  0).  Although  Ladanyi  et  al.  did  not  consider  ice  strengthening  due 
to  confinement,  the  shaded  zone  used  higher  ice  strengths  corresponding  to  ffc  =  10  MPa. 
The  MFS  upper  yield  stress  and  peak  strength  data  for  Dr  —  05%  are  from  Figs.  3.12  and 
3.18,  respectivdy. 

Given  the  uncertainties  in  the  ice  properties,  the  predicted  strengths  at  low 
confinement  compare  reasonably  well  with  the  measured  peak  strengths.  However,  there  is 
a  severe  problem  with  "strain  compatibility":  the  theory  requires  identical  strain  paths  for 
the  soil  skeleton  in  both  frozen  and  unfrozen  states,  and  the  pore  ice  cannot  be  broken. 
The  frozen  MFS  failed  at  Cp  =  7%  to  4%  (with  increasing  e),  whereas  unfrozen  sand 
requires  more  than  15%  strain  (test  22  in  Fig.  4.2)  and  unconfined  granular  ice  would 
probably  fail  at  less  than  1  or  2%. 

Comparison  at  high  confinement  is  less  satisfactory,  with  measured  peak  strengths 
being  much  larger  than  predicted  at  slow  and  moderate  strain  rates,  even  after  including 
higher  values  of  Qi  for  the  ice  matrix.  But  strain  compatibility  is  better  regarding  the 
sand  skeleton  since  frozen  MFS  failed  at  Cp  =  15%  to  23%  [Fig.  3.8(b)]. 

5.3  DISCUSSION 

Although  the  dilatancy— hardening  modd  represents  a  significant  advancement  in 
trying  to  understand  the  physical  mechanisms  controlling  the  behavior  of  dense  frozen 
sand,  the  results  in  Fig.  5.3  clearly  show  that  detailed  evaluation  of  the  model  requires 
experimental  data  covering  a  range  of  frozen  soil  conditions,  plus  data  on  undrained, 
unfrozen  sand  behavior.  Major  problems  with  the  modd  appear  to  be: 

1)  Strain  incompatibility  between  the  sand  skdeton  in  frozen  and  unfrozen  states. 

2)  The  assumption  that  the  sand  skdeton  has  identical  effective  stress  behavior  in 
frozen  and  unfrozen  states.  For  example,  the  forces  acting  on  sand  partides  will 
differ  significantly  when  surrounded  by  ice  rather  than  water  (see  bdow). 
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3)  Estimates  of  properties  for  pore  ice  having  grain  sizes  mnch  smaller  than  tested. 

4)  Although  Ladanyi  and  Morel  (1990)  imply  that  the  theory  applies  to 
"unconsolidated”  frozen  sand,  the  magnitude  of  the  effective  stress  initially  acting 
on  the  soil  skeleton  presumably  will  affect  frozen  soil  behavior.  For  example,  the 
data  in  Fig.  4.2  show  a  much  stiffer  response  for  dense  unfrozen  sand  with  increasing 
consolidation  (effective)  stress  even  though  all  samples  had  the  same  initial  density. 
Section  2.3.4  of  Andersen  (1991)  disciuses  item  2)  in  some  detail.  The  following 

summarizes  the  main  points.  The  mechanics  of  deformation  in  unfrozen  sands  that 
experience  no  particle  crushing  is  governed  by  the  magnitude  and  direction  of  the 
intergranular  contact  forces.  These  contact  forces,  when  integrated  over  a  large  area, 
describe  an  equilibrium  stress  state  in  terms  of  both  normal  and  shear  stresses.  In  contrast 
to  an  unfrozen  sand,  a  saturated  frozen  sand  is  a  continuous  medium  (i.e.,  the  pores 
contain  an  ice  matrix  that  has  a  significant  shear  strength).  This  means  that  both  the  sand 
particles  and  pore  ice  pairtidpate  actively  in  the  stress  transferring  process.  This  results  in 
significant  differences  between  the  types  of  forces  that  are  carried  by  the  sand  skeleton. 

In  a  frozen  sand,  the  portion  of  the  total  applied  stress  that  is  transmitted  by  the 
sand  skeleton  via  intergranular  contact  forces  may  not  form  an  equilibrium  stress  state  as 
is  the  case  in  unfrozen  sands.  This  is  illustrated  qualitatively  in  Fig.  5.4.  This  figure  is  a 
two-dimensional  free-body  diagram  of  an  individual  sand  partide  in  a  frozen  sand  matrix. 
The  partide  is  idealized  as  a  five-sided  polygon.  One  set  of  forces  acting  on  this  partide 
are  due  to  the  intergranular  contacts  with  adjacent  sand  grains  and  these  are  depicted  as 
solid  vectors.  Another  set  of  forces  acting  on  this  particle  are  those  transmitted  across  the 
ice-silicate  interface,  that  may  have  both  normal  and  shear  components  (Jellinek  1962). 
These  are  depicted  as  broken  arrows.  The  writers  recognize  the  inaccuracy  assodated  with 
this  schematic  diagram  in  the  assumption  that  the  five  intergranular  forces  are  coplanar. 
In  actuality,  there  are  probably  not  even  five  partides  in  contact  with  an  individual  grain 
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considering  all  thiee  dimensions.  This  presentation  is  only  given  to  discuss  the  idea  of 
particle  equilibrium  and  the  2— D  simplification  is  made  for  convenience. 

Since  ice  can  support  shear  stresses,  this  means  that  in  general,  neither  the 
intergranular  contact  forces  when  considered  separately,  nor  the  ice-silicate  interface  forces 
when  considered  separately,  will  form  an  equilibrium  force  condition  on  the  particle.  But 
when  these  are  considered  jointly,  and  if  the  particle  is  not  accderating,  the  net  resultant 
of  all  of  the  forces  acting  on  the  particle  must  be  zero  (neglecting  the  weight  of  the  sand 
particle),  i.e.,  the  total  forces  must  constitute  an  equilibrium  condition. 

Ice  will  creep  even  under  extremely  small  shear  stresses.  This  means  that  the 
magnitude  of  the  forces  transmitted  across  the  ice-silicate  interface  are  highly  time 
dependent  and  the  state  of  stress  in  the  pore  ice  will  eventually  tend  towards  hydrostatic. 
When  external  stresses  are  initially  applied  to  a  frozen  sand,  a  portion  of  that  stress  will  be 
carried  by  the  pore  ice.  Some  fraction  of  these  pore  ice  stresses  must  be  transmitted  across 
the  ice-silicate  interface.  Over  time,  if  the  strain  rate  goes  to  zero,  the  shear  stresses  being 
carried  by  the  ice  matrix  decrease  towards  zero  and  a  hydrostatic  state  of  stress  will  be 
approached  in  the  pore  ice.  For  this  limiting  condition,  both  the  ice-silicate  interface 
forces  (now  hydrostatic)  and  the  intergranular  contact  forces,  when  considered  separately, 
will  each  be  in  equilibrium  and  the  state  of  stress  determined  from  the  intergranular 
contact  forces  will  be  similar  to  that  in  an  unfrozen  sand  under  the  same  total  applied 
stresses.  Such  a  condition  may  only  apply  under  isotropic  stresses  after  consolidation  when 
the  strain  rate  is  very  close  to  zero.  In  general  however,  under  the  action  of  external 
stresses,  the  strain  rate  does  not  go  to  zero  and  there  will  always  be  some  shear  stresses 
acting  on  the  pore  ice  and  hence  the  intergranular  contact  forces  will  not  be  in  equilibrium 
when  considered  separately. 

If  the  intergranular  contact  forces  do  not  represent  a  state  of  equilibrium  at  the 
particle  level,  then  integrating  these  forces  over  a  large  area  will  generate  a  state  of  stress 
that  is  also  not  in  equilibrium.  Thus,  unless  the  intergranular  contact  forces  are  in 
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equilibrium,  the  state  of  stress  determined  from  these  will  not  be  similar  to  the  state  of 
stress  acting  on  the  same  sand  skeleton  in  an  unfrozen  state.  This  means  that  it  might  be 
possible  for  a  frozen  sand  to  undergo  a  strain  history  which  would  be  impossible  for  the 
same  system  in  an  unfrozen  state  (e.g.,  dilation  of  loose  frozen  sand  at  low  confinement  to 
less  than  Dr  =  0%).  In  the  writers’  opinion,  this  r^resents  a  significant  difference  between 
the  mechanics  of  particle  deformation  in  frozen  sands  as  opposed  to  unfrozen  sands. 


(«dw)  =  b 
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MIT  q-^  Plot  for  Ladanyi’s  Dilatancy  Hardening  Model 


state  Line  (SSL) 


Fignxe  5.2:  Critical  Confining  Pressure  Line  of  Unfrozen  MFS  from 

CIUC  Tests 


Figme  5.3:  Application  of  Ladanyi’s  Model  to  Dense  Frozen  MFS 

(Dr  =  95%,T  =  -9.5*C) 


FORC£  POLYGONS 


Sand  Intergranular  Contact  Forces 


Ice-Silicate  Interface  Forces 


Figure  5.4: 


Force  Equilibrium  of  Sand  Particle  In  Frozen  Sand 
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6.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  SUMMARY  AND  CONCLUSIONS 

Research  Objectives 

Development  of  rational  constitutive  idationships  to  modd  the  time-dependent 
strength-deformation  behavior  of  frozen  sand  for  use  in  design  practice  involves  three 
steps:  1)  a  complete  characterization  of  the  stress-strain— time  behavior  of  frozen  sand  as  a 
function  of  the  rdevant  experimental  variables  (i.e.,  sand  density,  confining  pressure,  strain 
rate  and  temperature);  2)  analysis  of  the  results  (often  using  deductive  reasoning  based  on 
comparisons  with  the  response  of  ice  and  unfrozen  sand)  to  identify  the  physical 
mechanisms  controlling  behavior;  and  3)  the  eventual  use  of  prindples  from  ice  mechanics, 
soil  mechanics,  composite  materials,  etc.  to  devdop  physically  based  constitutive  modds. 
This  research  focused  mainly  on  step  1)  because  an  extensive  literature  review  showed  that 
this  type  of  information  had  not  been  published  for  any  frozen  soil.  The  resulting 
experimental  data  set  represents  the  most  significant  contribution  of  the  research,  although 
progress  was  also  made  on  step  2). 

Experimental  Program 

A  high-pressure  low— temperature  triaxial  testing  system  was  developed  in  order  to 
measure  the  stress-strain  behavior  of  frozen  Manchester  Fine  Sand  (MFS)  from  very  small 
(0.01%)  to  very  large  (25%)  axial  strains.  This  system  incorporates  internal  force  and 
temperature  measurement,  lubricated  end  platens,  on-snedmen  axial  strain  measurement, 
pressure  control  to  10  MPa,  temperature  regulation  to  ^0.35*  C  and  volumetric  strain 
measurement  to  ±0.2%.  The  "prindpal”  testing  program  on  frozen  MFS  evaluated  the 
following  variables  for  spedmens  sheared  in  triaxial  compression  at  a  temperature  of 
T  = -9.5’ C  (Table  3.1): 

•  Relative  density  (Dr)  ranging  from  20  to  100%; 

•  Confining  pressures  (ctc)  of  0.1,  2,  5  and  10  MPa; 
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•  Axial  strain  rates  (c)  of  3  x  10‘*/sec  (slow),  3  x  lO'B/sec  (moderate)  and 
4  X  10’4/sec  (fast). 

A  "supplemental"  program  evaluated  the  influence  of  varying  the  temperature  from  to 
— 28*C  for  loose  sand  sheared  at  Cc  =  0.1  MPa  and  £  =  3  x  lO’Vsec  (moderate  rate). 

A  similar  high-pressure  triaxial  testing  system  was  developed  to  measure  the 
stress— strain  behavior  of  unfrozen  MFS  as  a  function  of  relative  density  and  effective 
(consolidation)  pressure.  This  program,  which  included  32  successful  consolidated— 
undrained  (CIUC)  and  6  successful  consolidated— drained  (CIDC)  triaxial  compression  tests 
(Table  4.1),  was  conducted  for  two  reasons:  1)  to  compare  the  strength-deformation 
behavior  of  the  same  sand  in  frozen  and  unfrozen  states;  and  2)  to  obtain  sand  properties 
required  for  application  of  Ladanyi’s  (1985)  dilatancy— hardening  model  that  attempts  to 
quantify  the  relative  contributions  of  the  frictional  resistance  of  the  soil  skeleton  and  of  the 
tensile  and  compressive  strength  of  the  ice  matrix. 

Most  of  the  specimens  used  for  both  the  frozen  and  unfrozen  triaxial  testing 
program,  were  prepared  by  multiple  sieve  pluviation  (MSP),  which  produced  very  uniform 
sand  specimens  having  a  vdde  range  of  relative  densities. 

General  Overview  of  Stress-Strain  Behavior 
of  Frozen  MFS 

The  use  of  on-specimen  axial  strain  measurements  provided  the  first  complete 
characterization  of  the  stress— strain  behavior  of  a  frozen  soil  from  very  small  to  very  large 
strains.  This  enabled  determination  of  the  following  stress-strain  parameters  from  plots  of 
axial  strain  (fa)  versus  deviator  stress  (Q  =  o-i  —  <73)  and  volumetric  strain  (ev  =  Av/vo), 
as  illustrated  in  Fig.  3.1. 

1)  Definition  of  small  strain  behavior,  such  as  values  of  Young’s  modulus  (E) 
representing  the  initial  slope  of  the  stress-strain  curve. 

2)  Values  of  the  upper  yield  stress  (Quy)>  which  represents  the  knee  of  the  stress-strain 
curve  where  very  significant  amounts  of  plastic  deformation  begin  to  occur. 
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3)  The  nature  of  the  post  upper  yield  behavior,  i.e.,  the  d^ree  of  strain  hardening  or 
strain  softening  after  the  knee.  Most  curves  exhibited  a  distinct  double  yield 
behavior  similar  to  that  of  fine-grained  polycrystalline  ice  within  the  ductile  to 
brittle  transitional  deformational  regime. 

4)  Values  of  the  peak  stress  (Qp)  and  post  peak  behavior. 

5)  The  volumetric  behavior,  including  the  rate  of  dilation  =  dfv/dea- 

The  general  shape  of  the  stress-strain  curves  was  classified  into  seven  different 
types.  These  are  shown  in  Fig.  3.2  which  plots  deviator  stress  divided  by  the  upper  yield 
stress  (Q/Quy)  for  tests  sheared  at  "low"  and  "high"  confining  pressures.  The  effects  of 
changes  in  relative  density,  confining  pressure,  strain  rate  and  temperature  on  stress-strain 
behavior  are  illustrated  in  Figs.  3.6,  3.7,  3.8  and  3.19. 

Small  Strain  Behavior  of  Frozen  MFS 

The  results  in  Fig.  3.3  show  that  the  initial  Young’s  modulus  (£)  equals  26.6  ±  4.5 
MPa  and  is  essentially  independent  of  changes  in  density,  confinement  and  strain  rate. 
Section  3.5  also  indicates  that  E  does  not  vary  with  temperature.  Section  3.3.2  shows  that 
the  initial  modulus  of  frozen  MFS  can  be  reasonably  predicted  by  treating  the  system  as  a 
composite  material  using  Counto’s  (1964)  model.  In  this  isostrain  model,  the  frozen  sand 
was  idealized  as  a  cube  of  "stiff"  sand  grains  embedded  in  a  "soft"  ice  matrix. 

Large  Strain  Behavior  of  Frozen  MFS 

The  magnitude  of  the  upper  yidd  stress  (Quy)  is  independent  of  sand  density  (Fig. 
3.10)  and  decreases  only  slightly  with  confining  pressure  (Fig.  3.11).  This  suggests  that 
Quy  is  not  affected  by  the  frictional  characteristics  of  the  sand  skeleton.  Although  Quy  is 
about  double  the  unconfined  peak  strength  of  fine-grained  polycrystalline  ice  at 
temperatures  near  — 10*C,  both  systems  exhibit  similar  sensitivities  to  strain  rate  (Fig. 
3.14).  This  observation  led  to  the  suggestion  in  Section  3.4.2  that  the  physical  mechanisms 
controlling  the  upper  yield  behavior  of  frozen  MFS  might  be  similar  to  those  controlling 
the  strength  of  ice.  However,  initial  results  from  low-density,  low-confinement,  moderate 
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strain  rate  tests  at  varying  temperatures  show  that  the  temperature  sensitivity  of  frozen 
MFS  is  very  different  from  that  predicted  for  ice  (Fig.  3.21).  Thus  the  above  hypothesis 
requires  further  experimental  evaluation,  such  as  proposed  in  Section  6.2. 

The  peak  strength  (Qp)  of  frozen  MFS  is  greatly  affected  by  the  frictional  resistance 
of  the  sand  skeleton,  except  for  Type  I  curves  wherein  failure  occurs  at  the  upper  yield 
stress.  This  Qp  =  Quy  condition  occurred  for  loose  sand  having  low  confinement  and 
sheared  either  at  a  fast  rate  or  at  a  low  temperature.  All  other  conditions  (lower  e  or 
higher  T,  higher  Dr,  and/or  higher  led  to  post  yield  strain  hardening  and  failure  at 
larger  strains.  Under  these  latter  conditions,  the  peak  strength  for  shearing  at  T  =  -9.5*  C 
was  shown  to:  1)  increase  linearly  with  relative  density,  with  the  rate  of  increase  becoming 
larger  with  confinement  (Fig.  3.15);  and  2)  increase  nonlinearly  with  confinement,  with  the 
rate  being  larger  at  higher  density  (Fig.  3.17).  These  strengthening  effects  are  less  rate 
sensitive  than  granular  ice  (Fig.  3.18)  and  presumably  reflect  a  greater  contribution  of  the 
frictional  resistance  of  the  sand  skeleton  to  the  peak  strength.  However,  initial  results  from 
tests  at  varying  temperature  indicate  that  the  temperature  sensitivity  of  Qp  for  frozen 
MFS  is  much  greater  than  that  for  granular  ice  at  temperatures  ranging  from  — 5*C  to 
-15*  C  (Fig.  3.21). 

Behavior  of  Unfrozen  MFS 

Section  4  presented  stress-strain  curves  and  effective  stress  paths  from  isotropically 
(hydrostatically)  consolidated— undrained(ClUC)  and  consolidated-drained  (CIDC)  triaxial 
compression  tests  that  illustrate  the  effects  of  changes  in  relative  density  and  confining 
pressure  on  the  behavior  of  unfrozen  MFS.  These  data  were  used  to  define  the  location  of 
the  large  strain  Steady  State  Line  (SSL)  in  void  ratio-log  stress  space  (Fig.  4.7)  and  then 
evaluated  in  terms  of  the  state  parameter  ^  (Fig.  4.6)  proposed  by  Been  and  Jefferies 
(1985).  This  approach  unified  the  combined  effects  of  changes  in  density  and  confining 
pressure  for  CIUC  tests  as  illustrated  by  plots  of  9  versus  the  pore  pressure  parameter  at 
failure  (Fig.  4.8)  and  undiained  strength  ratio  (Fig.  4.9). 
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Application  of  Dilatancy— Hardening  Model 
Section  5  summarized  the  dilatancy— hardening  model  developed  by  Ladanyi  (1985) 
to  predict  the  strength  of  frozen  dense  sand  in  terms  of  the  frictional  resistance  of  the  sand 
skeleton  and  of  the  tensile  and  compressive  strength  of  the  ice  matrix.  Application  of  this 
approach  to  frozen  dense  MFS  required  the  determination  of  the  critical  confining  pressure 
(^'scr)i  which  was  obtained  from  the  ^  versus  Af  relationship  at  Af  =  0  to  give  a'zcr  = 
2.37  MPa  (Fig.  5.2).  Figure  5.3  compared  predicted  and  measured  strengths  at  varying 
strain  rates.  At  low  confinement  (oc  =  0.1  MPa),  the  agreement  is  considered  reasonable 
given  the  uncertainties  in  the  assumed  ice  properties,  but  there  is  a  severe  problem  with 
"strain  compatibility"  as  discussed  in  Section  5.2.  At  high  confinement  (oc  =  10  MPa), 
the  measured  Qp  is  much  larger  than  predicted  at  slow  and  moderate  strain  rates.  Section 
5.3  discussed  other  potential  problems  with  the  dilatancy— hardening  model,  with  emphasis 
on  how  the  ice  matrix  might  afreet  the  stress-strain  behavior  of  the  sand  skeleton  which  is 
assumed  in  the  model  to  be  the  same  for  frozen  and  unfrozen  sand. 

6.2  RECOMMENDATIONS 

1)  The  results  in  Section  3.5  from  the  initial  study  of  the  efrects  of  changing 
temperature  on  the  behavior  of  frozen  MFS  indicate  trends  very  difrerent  from  those 
obtained  for  changing  strain  rate  (both  compared  to  their  efrects  on  ice  behavior). 
This  study  should  be  extended  to  include  a  fairly  comprehensive  program  of  testing 
over  a  range  of  experimental  conditions  such  as  outlined  in  Table  3.2. 

2)  The  undrained  (constant  volume)  stress-strain  data  for  unfrozen  dense  MFS  in 
Fig.  4.2  show  that  increases  in  the  preshear  effective  stress  cause  a  very  large 
increase  in  the  stiffness  of  the  sand  skeleton.  Such  changes  in  stiffness  presumably 
could  have  a  very  significant  affect  on  the  strength-deformation  behavior  of  frozen 
sand  at  strains  near  and  beyond  the  upper  yield  region.  Consequently,  tests  should 
be  conducted  on  sand  that  is  frozen  after  application  of  the  confining  pressure, 
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if 

rather  than  vice  versa  as  has  been  done  for  all  prior  experimental  programs  on 
firosen  soil.  This  will  require  a  more  complex  freezing  technique,  but  should 
simulate  more  closdy  actual  in  situ  conditions  for  "deep"  permafrost  and  for  soil 
stabilization  via  artificial  freezing.  Such  testing  also  might  help  to  further  evaluate 
Ladanyi’s  dilatancy— hardening  model. 

3)  No  information  exists  about  the  behavior  of  ice  having  very  small  grain  sizes  that 
presumably  exist  within  most  frozen  sands.  An  investigation  of  the  response  of 
frozen  spheres  having  large  changes  in  the  mean  particle  diameter  (e.g..  using  glass 
beads)  might  prove  useful.  Likewise,  changes  in  the  particle  surface  (roughness  and 
affinity  for  water)  may  shed  light  on  behavior  at  the  ice-silicate  interface. 

4)  The  modding  of  frozen  sand  as  a  composite  material  deserves  further  study  to 
determine  if  this  approach  can  be  extended  to  larger  strain  levels.  The  Discrete 
Element  Method  (DEC)  also  might  be  considered,  although  introduction  of  a 
"visoelastic"  matrix  will  certainly  make  such  analyses  highly  complex  compared  to 
current  capabilities. 
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APPENDIX  A:  TABULATED  DATA  FROM  TEST 
PROGRAM  ON  FROZEN  MFS 

Summaxy  of  Specimen  Physical  Properties  (2  pages) 

Measured  Strength  and  Deformation  Properties:  Principal  Program  (2  pages) 

Statistical  Analysis  of  Small  Strain  Data:  Principal  Program 

Summary  of  Small  Strain  Behavior:  Principal  Program 

Statistical  Analysis  of  Upper  Yield  Bdiavior:  Principal  Program 

Summary  of  Upper  Yield  Behavior:  Principal  Program 

Summary  of  Work  Hardening  Behavior:  Principal  Program 

Summary  of  Peak  Strength  Behavior:  Principal  Program  (2  pages) 

Summary  of  Volumetric  Behavior:  Principal  Program  (4  pages) 

Summary  of  Effect  of  Specimen  Preparation  Technique:  Principal  Program 
Effect  of  Temperature  on  Stress-Strain  Behavior 
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TABLE  A. 5  STATISTICAL  ANALYSIS  OF  UPPEE  YIELD  BEflAHOI 
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TABLE  A.6  SBHAKY  OF  UPPER  YIELD  BEHAVIOR 
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TABLE  A.7  SmOUIY  OF  VOU  HARDENING  BEHAVIOR 
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^  Premature  failure  before  reaching  "true”  peak. 

3  No  ACDT  measurements  therefore  values  are  approximate. 

5  Represents  end  of  test  measurements.  May  not  be  "true"  peak. 
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